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NRL  COMPACT  ACCELERATOR  THEORY  STUDIES 


INTRODUCTION 

This  report  contains  seven  short  papers  which  will  appear  in  the 
Proceedings  of  the  1989  Annual  DARPA/SDIO/Services  Charged  Particle  Beam 
Review  which  took  place  at  the  Naval  Postgraduate  School  in  Monterey,  CA 
during  18-21  September,  1989.  The  papers  describe  theoretical  studies  which 
support  two  compact  accelerator  experiments:  the  Spiral  Line  Induction 
Accelerator  (SLIA)  at  Pulse  Sciences,  Inc.  and  the  Recirculating  Linear 
Accelerator  (RLA)  at  Sandia  National  Laboratories. 

Both  accelerators  use  high  current  induction  accelerator  technology 
developed  for  the  ATA  and  RADLAC  accelerators.  To  reduce  the  weight  and 
volume  of  the  accelerator,  the  beam  will  be  recirculated  through  the  same 
accelerating  cavity  several  times.  The  SLIA  uses  ferrite-loaded  cavities 
similar  to  those  used  on  ATA  and  employs  stellarator  coils  to  provide  strong 
focusing  in  the  curved  sections  of  the  beamline.  The  RLA  uses  dielectric 
cavities  and  relies  on  the  electrostatic  focusing  of  an  ion  channel  to 
provide  strong  focusing.  The  NRL  studies  presented  here  include  analytical 
and  simulation  modeling  of  beam  transport  and  accelerator  instabili ties . 

The  simulations  have  employed  the  ELBA  code,  a  new  three-dimensional 
particle  code  capable  of  treating  both  straight  and  curved  sections  of  the 
accelerator. 

A  brief  summary  of  each  paper  and  a  list  of  co-authors  for  each  are 
provided  below. 

A.  ELBA  -  A  Three-Dimensional  Particle  Simulation  Code:  ELBA  is  a  beam 
simulation  code  which  treats  particle  dynamics  fully-relativistically  and 
solves  Maxwell's  equations  using  r,  9,  C  =  ct  -  z,  and  t  as  the  independent 
variables.  The  code  treats  propagation  in  curved  sections  of  an  accelerator 
by  using  a  mapping  technique  developed  by  Friedman  of  LLNL.  The  code  has 
been  used  to  model  beam  transport  and  three-wave  instbilitv  growth  in  the 
SLIA  and  inductive  erosion  and  off-axis  drifting  in  the  RLA.  (Joyce.  Krall, 
Slinker ) 

B.  Stability  Regimes  in  a  Rotating-Ouadrupole  Focusing  Accelerator: 

The  SLIA  will  use  rotating  quadrupole  or  orellerator  fields  to  provide 
strong  focusing  in  the  curved  sections  of  the  accelerator.  Beam  disruption 
from  the  three-wave  electromagnetic  instability  has  been  a  major  concern.  A 
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linear  dispersion  relation  has  been  derived  for  this  instability  assuming  a 
thin,  cold  beam  propagating  in  a  cylindrical  waveguide  with  specified 
solenoidal  and  rotating  quadrupole  external  magnetic  fields,  The  beam 
centroid  may  be  orbit  unstable,  three-wave  unstable  or  fully  stable 
depending  on  the  beam  and  external  field  parameters.  Analytical  expressions 
for  the  boundaries  between  these  regimes  have  also  been  derived,  making  it 
possible  to  illustrate  graphically  the  regions  of  stability  and  instability. 
The  model  predicts  that  the  SLIA  proof-of-concept  experiment  should  be  able 
to  avoid  the  instability  entirely  if  operated  in  the  reversed-helici ty 
regime  suggested  by  Chernin  and  Hughes.  At  higher  energies,  it  should  be 
possible  to  remain  in  the  stable  regime  by  decreasing  the  quadrupole 
gradient  in  the  high  energy  bends.  These  results  are  consistent  with 
independent  calculations  by  MRC  and  SAIC.  (Tang,  Sprangle,  Krall,  Serafim, 
Mako) 

C.  Three-Dimensional  Simulations  of  Electromagnetically  Stable  and 
Unstable  Regimes  in  the  Spiral  Line  Induction  Accelerator:  The  ELBA 
simulation  code  has  been  used  to  simulate  the  Three-wave  instability 
described  above.  The  simulations  included  the  appropriate  external  fields 
and  employed  the  STELMAT  algorithm  from  SAIC  for  matching  the  beam 
initially.  Five  simulations  were  carried  out  with  y  =  7,  7  kA  beam  current, 
200  G/cm  quadrupole  gradient,  0.5  cnT^  quadrupole  wavenumber,  3  cm  wall 
radius,  0.158  cm-rad  normalized  rms  emittance,  and  a  guide  field  from  1 
to  6  kG.  The  linear  theory  described  above  predicts  three-wave  instability 
for  Bz  =  1,  2  and  4.5  kG,  orbit  instability  at  B^  =  4  kG,  and  stability  at 
Bz  =  9  kG,  and  the  individual  simulations  confirm  those  predictions.  The 
simulation  growth  rates  in  the  three-wave  unstable  regime  were  significantly 
lower  than  those  predicted  by  the  dispersion  relation,  however,  and  the 
source  of  this  discrepancy  is  still  under  investigation.  Stable  propagation 
was  also  observed  for  Bz  =  -5  kG  and  a  low  quadrupole  gradient,  but  this 
stability  regime  is  very  narrow  and  thus  much  less  useful  experimentally. 

The  results  support  the  conclusion  from  the  analytical  models  that  the  SLIA 
proof-of-concept  experiment  can  operate  entirely  in  the  stable  regime. 
(Krall,  Tang,  Joyce) 

D.  Three-Dimensional  Simulations  of  Beam  Transport  and  Emittance  Growth 
in  the  Spiral  Line  Induction  Accelerator:  The  ELBA  code  was  used  to  model 
transport  in  an  idealized  transition  region  between  purely-solenoidal  and 
solenoidal/stellarator  sections  of  SLIA.  Emittance  growth  in  these 
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transitions  is  a  major  potential  problem  in  some  parameter  regimes.  When  a 
high-eccentricity  equilibrium  was  chosen  in  conjunction  with  a  long 
stellarator  pitch  length,  normalized  rms  emittance  grew  by  as  much  as  47%  in 
a  single  transition.  When  the  pitch  length  was  shortened  to  a  length 
comparable  with  the  ramp  length  of  the  stellarator  field,  emittance  growth 
was  reduced  substantially.  The  choice  of  a  low  eccentricity  equilibrium 
resulted  in  negligible  emittance  growth,  but  this  regime  may  be  less 
desirable  for  other  reasons.  (Krall,  Joyce) 

E.  NRL  Studies  of  RLA  Accelerator  Physics:  The  Recirculating  Linear 
Accelerator  (RLA)  currently  under  development  at  Sandia  National 
Laboratories  uses  ion-focused  regime  transport  (IFR)  to  guide  the  beam.  In 
a  racetrack  configuration,  the  ion  channel  must  survive  for  as  long  as  400 
nsec,  a  time  long  compared  with  the  characteristic  time  for  a  bare  ion 
channel  to  expand  due  to  its  own  space  charge.  The  FRIEZR  code  has  been 
modified  to  treat  long-time-scale  ion  motion.  PLasma  electron  and  ion 
dynamics  are  included,  and  an  artificially  low  ion  mass  (hydrogen  instead  of 
xenon)  is  used  with  an  appropriately  rescaled  pulse  length.  Je  first 
investigated  whether  plasma  electrons  from  the  wall  could  neutralize  the  ion 
space  charge  and  slow  the  expansion  of  the  channel  after  the  beam  pulse  has 
ended;  however,  the  inclusion  of  these  electrons  was  not  very  effective. 
Somewhat  more  promising  results  were  obtained  when  the  beam  created  part  of 
the  channel  itself,  thus  replenishing  part  of  the  channel  after  each 
recirculation.  This  strategy  could  be  employed  by  raising  the  background 
gas  pressure  to  a  few  mtorr.  We  also  made  analytical  estimates  of  the 
characteristic  time  for  a  bare  channel  to  drift  off-axis  due  to 
electrostatic  image  charges  in  the  wall.  Although  this  effect  is  small,  it 
cannot  be  readily  corrected  for  since  the  electrostatic  equilibrium  in 
inherently  unstable.  (Joyce,  Hubbard,  Fernsler) 

F.  Erosion  of  Electron  Beams  in  the  RLA  Accelerator:  The  ELBA 
simulation  code  has  been  used  to  study  inductive  erosion,  beam  drift  and 
current  loss  in  both  straight  and  curved  sections  of  the  RLA.  The  observed 
erosion  rate  in  straight  sections  agrees  well  with  analytical  estimates  and 
2-D  FRIEZR  simulations.  In  curved  sections,  the  erosion  rate  is  typically 
30-50%  higher.  If  the  vertical  field  By  is  well-matched  to  the  nominal  beam 
energy  and  bend  radius  of  curvature,  inductile  energy  loss  at  the  beam  head 
may  causes  the  beam  centroid  in  that  region  to  drift  to  the  inside.  Plasma 
electrons  are  ejected  along  the  vertical  field  lines,  but  this  complicated 
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process  does  not  appear  to  degrade  the  beam  significantly.  A  10 X  mismatch 
in  the  vertical  field  causes  modest  drifts  in  the  beam  centroid  and  has  a 
small  effect  on  the  erosion  rate. 

Beam  transport  simulations  of  a  ’’recirculated"  RLA  beam  have  been 
carried  out  by  using  an  expanded  and  weaker  channel  and  a  higher  beam 
energy.  As  expected,  the  inductive  erosion  rates  are  small,  but  significant 
drifting  in  the  beam  centroid  is  often  observed  because  the  channel  does  not 
provide  a  strong  centering  force.  (Hubbard,  Joyce,  Slinker) 

G.  Atomic  and  Plasma  Surface  Interactions  in  Accelerators  and  IFR 
Cells:  This  paper  summarizes  the  various  atomic  physics  processes  which  may 
occur  in  accelerators  and  IFR  cells.  Beam  impact  ionization  of  the  ambient 
or  outgassed  species  is  negligible  for  the  SLIA  but  may  be  important  for  the 
betatron  and  RLA.  Accelerated  plasma  electrons  and  ions  may  impact  the 
wall,  creating  a  substantial  surface  plasma.  In  addition,  fast  neutral 
particles,  creating  by  charge  exchange  with  plasma  ions,  may  also  impact  the 
wall  surface.  In  induction  accelerators,  these  plasmas  can  lead  to 
undesirable  electrical  breakdown  in  the  accelerating  gaps,  a  process  which 
has  been  observed  experimental  in  the  laser-guided  ATA  accelerator.  Both 
analytical  and  ab-initio  computer  models  are  currently  being  explored  to 
treat  these  processes.  (Ali) 
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ELBA  -  A  Three-Dimensional  Particle  Simulation  Code* 


Glenn  Joyce,  Jonathan  Krall,  Steve  Slinker 

Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory 


INTRODUCTION: 

The  study  of  the  interaction  of  relativistic  electron  beams  with 
plasmas  and  with  external  magnetic  field  configurations  has  required  the 
development  of  sophisticated  particle  simulation  codes.  We  developed  the 
FRIEZR  code  to  look,  at  axisymmetric  problems  of  beam  propagation  in  the  Ion 
Focus  Regime  (IFR).  A  number  of  investigations  were  made  using  the  models 
integrated  into  that  code.  Ve  are  new  progressing  to  the  study  of  three- 
dimensional  phenomena  associated  with  beam  propagation  through  diffuse 
plasmas  and  with  beam  transport  in  compact  accelerators.  We  have  developed 
a  three-dimensional  code,  ELBA,  to  investigate  the  physics  of  some  of  these 
phenomena.  The  ELBA  code  is  based  largely  on  FRIEZR  with  the  addition  of 
three-dimensional  electromagnetic  fields.  It  is  formulated  in  the 
(r,9,C=ct-z)  coordinate  system,  has  fully  relativistic  beam  and  plasma 
particles,  and  is  fully  electromagnetic.  In  addition,  we  have  included  the 
capability  of  simulating  the  transport  of  particle  beams  around  bends  in 
the  limit  of  a  large  radius  of  curvature.  The  latter  formulation  is  based 
on  a  mapping  scheme  developed  by  Alex  Friedman.  The  mapping  algorithm  does 
not  depend  on  the  radius  of  curvature,  but  the  electromagnetic  fields 
include  an  expansion  to  first  order  in  the  inverse  aspect  ratio  of  the  the 
simulation  waveguide  radius  to  the  radius  of  curvature. 

GENERAL  ISSUES: 

Because  of  the  nature  of  high  energy  electrons,  a  number  of  conventions 
have  been  adopted  for  their  study  which  are  not  usually  found  in  particle 
simulation  codes.  It  is  convenient  to  replace  the  space  and  time  variables 
z  and  t  by  C  =  ct-z  and  t.  Relativistic  beam  particles,  in  terms  of  these 
variables,  exhibit  very  slow  motion  in  the  axial,  C,  direction.  A  beam 
particle  moving  exactly  at  the  speed  of  light  would  have  a  constant  value 
of  C.  The  actual  motion  will  result  in  a  slow  increase  in  C  position  of 
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the  particle  depending  on  its  energy  and  perpendicular  velocity, 
plasma  particles  on  the  other  hand  appear  to  be  streaming  by  the  beam  at 
roughly  the  speed  of  light.  Ve  are  usually  interested  primarily  in  the 
beam  dynamics  so  the  C,  t  coordinates  are  generally  veil-suited  for  our  use. 
A  second  consequence  of  the  C.  t  coordinates  concerns  the  boundary 
conditions  for  the  electromagnetic  fields.  No  information  can  flov  forvard 
in  these  variables  so  the  propagation  of  all  electromagnetic  vaves  is 
toward  larger  C-  As  a  result,  ve  can  predetermine  the  values  of  the 
electromagnetic  fields  at  the  front  end  of  the  simulation  system  and  they 
will  remain  unchanged  by  the  beam. 

The  fact  that  ve  do  not  have  to  treat  information  moving  in  the  forvard 
direction  allows  us  to  handle  particle  data  in  a  particularly  simple  way. 
The  simulation  region  is  divided  axially  into  a  number  of  "slices"  with  the 
slice  length  corresponding  to  the  grid  sice  along  the  C  axis.  The  point 
C=0  is  generally  associated  with  initial  position  of  the  beam  head  and 
C=C_  „  is  associated  with  the  beam  tail.  At  any  given  time  each  particle, 
whether  it  is  a  beam  or  plasma  particle,  is  associated  with  a  particular 
slice.  As  time  progresses,  the  particle  either  stays  in  that  slice  or 
moves  to  the  next  slice  at  larger  C-  A  particle  at  any  given  slice  cannot 
be  influenced  by  a  particle  in  another  slice  which  is  nearer  the  tail  beam. 
As  a  result,  not  all  particles  need  to  be  in  the  computer  simultaneously. 

At  each  time  step  ve  can  begin  with  tne  first  slice,  read  the  particle 
variables  for  this  slice  (with  some  associated  field,  charge,  and  current 
density  quantities),  calculate  the  electromagnetic  fields  for  the  slice, 
push  the  particles,  calculate  the  charge  and  current  densities,  store  the 
particle, density,  and  field  information  on  disk  and  begin  with  the  next 
slice  working  our  way  from  the  front  to  the  back  of  ’the  beam.  As  particles 
move  from  slice  to  slice,  they  are  marked  as  having  moved  to  next  slice  and 
stored  appropriately.  The  marking  process  requires  that  we  must  do  a 
particle  sort  at  each  time  step  anc  at  each  slice.  The  sort  is  into  two 
regions:  the  particles  which  remain  in  the  slice,  and  those  which  move  to 
the  next  slice.  At  the  time  that  ve  do  the  sort,  we  also  handle  particles 
which  have  crossed  the  radial  simulation  boundary.  Normally,  ve  excise 
these  particles  from  the  simulation.  As  a  result  of  this  buffering  system, 
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we  can  simulate  beams  wi th  an  almost  unlimited  number  of  particles  although 
ve  pay  a  price  in  I/O.  A  single  rime  step  represents  the  processing  of  all 
the  slices  in  the  simulation  region  from  the  beam  head  to  the  beam  tail. 

No  particle  can  move  to  a  slice  nearer  the  beam  head  since  the  slice 
"velocity"  is  c.  The  beam  particles  slovly  move  from  slice  to  slice  while 
the  plasma  particles  change  slices  frequently.  In  order  to  make  the 
bookkeeping  simple,  ve  choose  the  time  step  to  be  smaller  than  uC-'Zc.  This 
constraint  is  also  imposed  on  us  by  the  field  solver. 

ELECTROMAGNETIC  FIELDS : 

Ve  solve  the  electromagnetic  fields  from  a  representation  of  the  scalar 
and  vector  potentials.  The  equations  are  solved  in  muc.n  the  same  manner  as 
FRIEZR  using  the  gauge  7,*A,*0.  Here.  A.  is  the  component  of  the  vector 
potential  perpendicular  to  the  C  axis.  The  fields  are  solved  in  much  t.ne 
same  vay  as  in  FRIEZR.  Ve  solve  mode  by  mode  in  a  Fourier  transformed 
space  vhere  the  Fourier  transform  is  taken  in  the  azimuthal  direction.  Ve 

vill  discuss  the  electromagnetic  field  solution  in  more  detail  in  a  more 
complete  paper  on  the  code. 

Ve  have  chosen  a  radial  grid  which  is  nonlinear  in  the  r  coordinate. 

Since  the  beam  particles  normally  lie  near  the  axis  of  t.ne  system,  i:  is 

convenient  to  have  the  grids  more  closely  spaced  here  than  near  the  veils. 

An  evenly  spaced  grid  in  the  coordinate  u  =  in(l  -  r“)  has  that  property. 

2 

Near  the  axis  u  -  i  ,  which  approximates  a  constant  current  density  grid, 
while  for  large  r  the  grid  spacing  increases  exponentially.  Our  experience 
is  that  this  grid  works  quite  well,  although  we  must  be  careful  in  its 
implementation. 

INITIAL  CONDITIONS: 

Although  the  choice  of  the  (C, t)  coordinate  system  simplifies  the 
boundary  conditions  at  C=0,  some  decisions  must  be  made  about  the 
laboratory  z=0  and  the  initial  conditions.  The  simplest  assumption  is  that 
it  is  not  important  where  the  beam  is  initially.  In  this  case,  the  beam 
cnarge  can  be  turned  on  adiabatically,  and  the  long  term  behavior  of  the 
system  can  be  studied.  Usually,  however,  we  want  to  simulate  the  beam 
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entering  a  regi~r,  of  interest  from  some  other  region  and  ve  must  decide  ho 
to  treat  (z--0.t-0).  After  a  number  of  trials,  ve  have  i  mp  1  emen :  eb  the 
following  initialisation  steps.  At  t=0.  ve  assume  that  the  head  of  the 
be'  r  is  at  s=0  and  that  the  rest  of  the  beam  is  out  of  the  region  of 
interest.  For  initial  conditions,  ve  choose  a  time  independent  solution 
for  the  electromagnetic  fields  appropriate  to  a  beam  in  a  vacuum.  The 
plasma  is  chosen  to  be  charge  neutral  and  cold.  As  the  beam  enters  the 
simulation  region  (z>0),  ve  push  the  beam  particles,  but  for  the  portion  o 
the  beam  remaining  outside  the  simulation  region,  ve  freeze  motion  in  the 
radial  direction  and  allow  the  beam  to  free  stream  in  the  C  direction.  In 
order  to  avoid  discontinui ties  at  the  z=0  boundary,  ve  allow  the  plasma  to 
react  to  the  fields  in  both  regions  of  the  simulation.  This  scenario 
provides  a  reasonable  model  for  the  entrance  of  the  beam  into  a  transport 
sys  tern. 

CONCLUSIONS: 

The  ELBA  code  has  proved  to  be  a  useful  tool  in  a  number  of  studies  of 
relativistic  beam  transport  and  propagation.  The  coordinate  system  limits 
its  applicability  to  beams  which  are  moving  at  roughly  the  speed  of  light. 
Ve  have  restricted  the  boundary  conditions  to  be  those  of  a  perfectly 
conducting  cylinder  so  that  the  electromagnetic  fields  can  be  solved 
quickly.  Ve  have  usually  opted  for  speed  ana  simplicity  over  generality. 
The  inclusion  of  bends  is  still  in  an  experimental  stage  in  that  ve  are 
continuing  to  study  and  make  improvements  in  the  algorithm. 

*Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA  Orde 
No.  4395,  Amendment  SO,  monitored  by  Naval  Surface  Warfare  Center. 
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Stability  Regimes  in  a  Rotating  Qnadrupole  Focusing  Accelerator' 

C.  M.  Tang,  P.  Sprangle,  J.  Krall,  P.  Serafim*  and  F.  Mako* 

Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  DC  20375-5000 

A  number  of  recent,  high  current  accelerator  configurations  utilize  strong  focusing 
Helds.  These  fields,  consisting  of  a  rotating  quadrupole  field  (or  stellarator  field)  and 
an  axial  guide  field,  increase  considerably  the  energy  mismatch  tolerance  of  the  device 
and  provide  confining  forces  against  the  beam  space  charge  forces.1,2  Two  such  devices 
are  the  modified  betatron  accelerator3  and  the  spiral  line  induction  accelerator  (SLIA).4 
The  SLIA,  to  be  constructed  bv  Pulse  Sciences,  Inc.,  will  utilize  the  strong  focusing  for 
transport,  along  the  curved  sections  of  the  beam  line  between  the  accelerating  cavities. 

The  use  of  strong  focusing  fields  has  a  potential  difficulty  in  that  they  can  lead  to 
various  types  of  beam  instabilities.5-'  We  find  that  the  electron  beam  centroid  can  be  i) 
orbit,  unstable  independent  of  the  waveguide  modes,  ii)  three-wave  unstable  or  iii)  fully 
stable.  The  conditions  for  each  regime  are  presented  in  the  limit  of  zero  beam  current. 
Algebraic  expressions  for  the  growth  rate  in  each  of  the  three- wave  unstable  regimes  are 
obtained,  and  they  are  in  general  agreement  with  the  expressions  given  in  Ref.  8. 

DISPERSION  RELATION.  The  external  magnetic  field  consists  of  an  axial  guide  field 
B0 .  and  a  rotat  ing  quadrupole  field  ( Bqz ,  B,VJ  ),  where 

BqT  =  —  Bqkq(x  sin  kqz  —  y  cos  kq: ),  Bqy  -  Bqkq\ x  cos  kq  z  +  y  sin  kqz ).  lie  -  />) 

Bqkq  is  the  quadrupole  gradient,  and  kq  is  the  wave  number  of  the  quadrupole  field. 
The  representation  for  the  quadrupole  field  in  Eqs.  (la-b)  is  valid  near  the  z-axis.  i.e.. 
(.r 3  +  ./2)1/2  v<  \q/2n.  In  equilibrium,  the  electron  beam  travels  along  the  axis  of  a 
circular  waveguide  at  velocity  v0  and  is  monoenereetic  with  7^  =  ( 1  —  ) ~ 1  ,  where 

30  =  r„/c. 

We  assume  the  electron  beam  propagates  within  a  perfectly  conducting  cylindrical 
waveguide  of  radius  rg .  We  include  induced  image  charges  and  currents  due  to  the  displaced 
beam.  Both  the  beam  radius  and  beam  centroid  displacement  are  assumed  to  be  small  in 
comparison  to  the  waveguide  radius. 

We  expect  TEu  mode  to  have  the  largest  growth  rate,  because  its  electric  field  peaks 
on  axis.  Its  vector  potential  can  be  written  as  A  =  An(r.tf)exp{i(fcr  —  ud)}  +  c.c.,  where 
uJ  is  the  radian  frequency,  k  is  the  wave  number  and  c.c.  denotes  the  complex  conjugate. 
The  boundary  condition,  requiring  the  tangential  component  of  the  electric  field  to  vanish 
on  the  perfectly  conducting  waveguide  surface,  gives  the  condition  J[(n\\rq)  =  0.  where 
/» 1 1  rq  is  the  smallest  positive  zero  of  Bessel  function  J[. 

Work  supported  by  Defense  Advanced  Research  Projects  Agency,  ARPA  Order  No.  -1305. 
Amendment  80.  monitored  by  Naval  Surface  Warfare  Center. 

1  Northeastern  University.  Boston,  MA  02115 
;  EM  Technologies.  Inc.,  Alexandria,  VA  22304 
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When  the  wave  equation  is  solved  simultaneously  with  the  particle  dynamics  equation, 
we  obtain  the  following  dispersion  relation 

WrWiWuWa  =  k2b  |(A-  +  kq  -  u >/v0Y  D_Wr  +  (k  -  u>/v„)~  D+  IF,  ,  (2) 

where  A-g  =  2(/{,/l7)/.tj1/(70(//21  r2  -  1  )J2(H\  i  rg )),  Ib  is  the  beam  current  in  kA,  U'r  - 
<m2/c2  —  k2  —  pj,  and  H/  -  u>2 / c2  —  (k  +  kq)2  —  n\x  are  contributions  from  the  right  - 
hand  circularly  polarized  (RHCP)  and  left-hand  circularly  polarized  (LHCP)  TE u  wave¬ 
guide  modes  respectively,  W,  =  K2  -  {dj  +  d2)  accounts  for  the  two  stable  beam  modes. 
H,ru  —  K2  —  { d \  —  d\)  admits  the  two  potentially  unstable  beam  modes,  d2r  —  A?  +  A 2. 

d\  =  ((A*!  4-  A'r/2)2  -  (A',4  -  A’4)) 1/2  ,  I<  =  k  -  +  kq/2.  D±  =  A'2  t  A,  A  -  A  A 

Ai  =  A0  -  A-g,  A |  =  (A0  -  kq/2)kq/2  -  A-2,  A'|  =  A7A-,,  A'*  =  |e|50//i„'>„m0r2  is  the 
relativistic  cyclotron  wave  number  associated  with  the  axial  held.  A",,  =  |»  |B7//i0-)0m„c2  is 
the  relativistic  cyclotron  wave  number  associated  with  the  rotating  quadrupole  hehl  and 
k,  =  (2(/i,/17)/(^272r2  ))A2 .  Equation  (2)  agrees  with  the  dispersion  relation  in  Ref.  (»i) 
with  the  vertical  field  set  to  zero. 

STABILITY  REGIMES  AND  ANALYTICAL  GROWTH  RATE  EXPRESSIONS.  The 
dispersion  relation.  Eq.  (2),  contains  i)  a  region  of  orbital  instability,  ii)  two  regions  of 
three- wave  instability  and  iii)  two  regions  of  stability.  The  stability  diagrams  are  obtained 
in  (kq,K0)  space  for  given  values  of  7.  rg  and  Bqk,r 

i)  Orbit  Lhistable  Regime.  The  electron  beam  in  this  configuration  can  be  orbit  and 
three-wave  unstable  when  (d\  -  d2)  _  0.  The  unstable  values  of  A„  are 

Acrj/,2  —  A.'  q  j  2  —  2I\q  A„  ■_  ArrO.3  =  kq;  2  -I-  2  A  q  .  ( •)  ) 

where  I\.q  —  \e\Bq/30^0m0c2  and  K0  =  |ed Bn,' 30-)0m0c2 .  Equation  (3)  is  in  agreement 
with  the  stability  condition  of  Ref.  1  in  the  limit  of  st  raight  cylindrical  geomet  ry  and 
zero  beam  current,  the  stability  condition  of  Ref.  2  in  the  limit  of  perfectly  conducting 
walls,  and  the  condition  for  beam  envelope  stability  of  Ref.  9  in  the  limit  of  zero  space 
charge. 

h)  Three- Wave  Unstable  Regimes.  Numerical  solutions  of  the  dispersion  relation  (2) 
indicate  that  the  three- wave  instability  occurs  when  the  RHCP  waveguide  mode  [or 
the  LIIC'P  waveguide  mode]  intersects,  in  the  (tu.A:)  plane,  the  appropriate  beam 
mode  given  bv  IT„  =  0  and  (d\  —  d2)  0.  For  A„  A’,;rif the  three- wave  is 

unstable  (Region  I)  when  the  RHC'P  waveguide  mode  intersects  the  beam  line  uy  r„  = 
(A-  +kq/2)+  \/d]  —  d\.  For  I\0  ■  Acr>/,3.  the  three-wave  is  unstable  (Region  II)  when 
the  RHCP  waveguide  mode  intersects  the  beam  line  u :.v0  =  (A’  -f  kq,  2)  -  \  d2  -  d2. 

We  obtain  algebraic  expressions  for  the  peak  growth  rates  /n?(^u.’/c).  where 
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where  u>\  =  y/k2  +  /xj, c,  u>2  =  t’0(A:  -1-  kq/2),  Au>2  =  va y/d]  -  d\,  Z>+U= = 
(Au^/vo)2  ±  A’i( Au)2/i’0)  —  A',  and  the  top  and  bottom  signs  in  (4)  refer  to  the 
three-wave  unstable  regions  I  and  II  respectively. 

iii)  Three-Wave  Stable  Regime  for  K„  <  7Cr-,»  -> ■  Stability  is  achieved  when  the  wave¬ 
guide  cut  off  frequency  fin c  is  sufficiently  large  so  that  intersection  with  either  of 
the  beam  lines,  defined  by  U"u  =  0  cannot  be  achieved.  The  three-wave  is  stable  if 

<7Pn  >  kq  +  2  {d\  -  d\ )  "  ,  where  q  =  (4 /(y2  -  2))  "  .  For  kq  <7pn,  and  defining 

(  =  k~(  1  +8K2/k2f)  —  2/  and  /  =  qptn(kq  -  g^u/2),  the  electron  beam  is  stable  for 
the  following  situations: 

a)  for  /  >  0  and  (  >  0,  the  stable  range  of  Ka  is  given  bv 

Acni.i  —  kq/2  y  "/  2  ■  A  r>  foe) 

b)  for  /  <  0  and  >  0,  the  stable  values  of  K0  are 

K0  <:  A'cr!(,i  =  smaller  of  (  kqj  2  -  C1 ,2/2.  K,-r,t.2  )  •  (  5M 

c)  for  /  <  0  and  4  <  0,  all  values  of  Ka  Afr,,,2  are  stable. 

iv)  Three- Wave  Stable  Regime  for  A~„  >  h^^}i  i .  The  three-wave  interaction  is  also  sta¬ 
ble  when  the  RHCP  waveguide  mode  does  not  intersect  ic  r„  —  (k  - f-  kqj  2)  —  \  d\  - 

for  A’o  >  Kent,  3-  This  occurs  when  _  A'7  -  2  (</j  -  d\)  ”  .  for  which  there  are 

two  stable  cases: 

a)  for  kq  <  <7/tn,  the  stable  range  of  A„  is  A\,  A'cro.3> 

b)  for  A:,  >  <7/zu,  the  stable  range  of  A’„  is 

A„  '•  I\Cnt,i  =  kq/ 2  -+-  4'1  ‘  / 2.  (6) 

In  the  limit  of  small  quadrupole  gradients,  large  yn  and  A*c  '  AUru.3-  the  stability 
condition10  is  approximately  K0  kq  — 

The  various  operating  regimes  are  illustrated  as  functions  of  kq  and  A’0  in  Fig.  1. 
for  7„  =  5,  rg  =  3  cm  and  quadrupole  gradient  Bqkq  =  200  G/cm.  Since  the  stability 
boundaries  are  obtained  in  the  limit  of  zero  beam  current,  the  area  of  the  two  stable 
regions  will  shrink  slightly  as  the  current  is  increased. 

To  operate  in  the  stable  regime  with  B„  0.  both  the  maximum  allowable  quadrupole 
gradient  and  the  quadrupole  wave  number  must  decrease  for  increasing  beam  energy.  Fig¬ 
ure  2  is  a  plot  of  maximum  quadrupole  gradient  versus  for  a  fixed  guide  field  B„  —  5  kG, 
rg  =  3  cm  and  two  different,  values  of  quadrupole  wave  number:  kq  =  0.1  cm-1  ( — )  and 
kq  =  0.05  cm  -1(-  -  - ). 

CONCLUSIONS.  Analytic  stability  boundaries  and  algebraic  expressions  of  the  growth 
rate  are  in  agreement  with  numerical  solutions  of  the  dispersion  relation.  These  results 
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suggest  that  the  SLIA  proof  of  concept  experiment  can  operate  entirely  in  the  stable 
regime. 


Fig.  1.  Plot  of  the  various  operating  Fig.  2.  Plot  of  maximum  quadrupole 

regimes  for  =  5,  rg  =  3  cm  and  quadrupole  gradient  versus  j0  for  guide  field  Ba  = 

gradient  Bqkq  =  200  G/cm.  5  kG  and  rg  =  3  cm. 
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THREE-DIMENSIONAL  SIMULATIONS  OF  ELECTROMAGNETICALLY  STABLE  AND 


UNSTABLE  REGIMES  IN  THE  SPIRAL  LINE  INDUCTION  ACCELERATOR* 


J.  KRALL,  C.  M.  TANG  and  G.  JOYCE 

Naval  Research  Laboratory 
Beam  Physics  Branch 
Plasma  Physics  Division 
Washington,  DC  20375-5000 

Magnetic  transport  of  high-current  electron  beams  is  an  important 
problem  in  induction  accelerator  physics.  One  such  accelerator  is  the 
spiral  line  induction  accelerator  (SLIA)*,  which  uses  solenoidal  focusing  in 
the  straight  sections  and  both  solenoidal  and  stellarator  (helical 
quadrupole)  focusing  in  the  curved  sections.  The  stellarator  fields,  which 
focus  the  space  charge  of  the  beam  over  a  large  energy  bandwidth,  make  the 
high-current  beam  subject  to  electromagnetic  instabilities,  such  as  the 
three-wave  instability. 

For  the  present  numerical  study  we  use  the  ELBA  code,  a  three- 
dimensional  particle  code  which  simulates  a  beam  propagating  within  a 
cylindrical  metallic  pipe.  The  full  set  of  Maxwell's  equations  along  with 
the  full  relativistic  motion  of  the  beam  particles  are  included.  The  beam 
is  initially  matched  to  the  field  configuration  to  minimize  initial 
oscillations.  This  matched  configuration  is  calculated  by  STELMAT3,  which 
was  adapted  for  use  in  ELBA. 

ELECTROMAGNETIC  THREE-WAVE  INSTABILITY 

The  electromagnetic  three-wave  instability  has  been  studied  by  various 
authors. Of  interest  here  are  theoretical  predictions  of  three-wave 
stable  field  configurations. In  the  limit  of  low  current  ( I ^  -»  0), 
analytic  expressions  for  stability  boundaries  in  (kq,Bz)  space  have  been 
obtained^,  where  kq  =  4rt/Xp  is  the  wavenumber  of  the  helical  quadrupole 
fields,  Xp  is  the  quadrupole  pitch  length,  and  Bz  is  the  axial  magnetic 
field.  These  boundaries  are  functions  only  of  the  amplitude  of  the 
quadrupole  gradient,  Bqkq,  the  relativistic  factor  for  the  beam,  y,  and  the 
waveguide  (pipe)  radius,  rg.  Stability  boundaries  are  plotted  in  Fig.  1  for 
parameters  of  interest  for  the  proposed  SLIA  proof-of-concept  experiment 
(PoCE):  Bqk„  =  200  g/cm,  y  =  7,  and  rg  =  3  cm.  The  operating  point  for  the 
PcCE  is  indicated  by  an  "x"  in  the  figure.  Note  that  the  diagram 
encompasses  two  three-wave  unstable  regimes,  which  are  physically  distinct, 
an  orbit  unstable  regime  and  two  stable  regimes.  Note  that  the  condition 
for  orbit  instability  in  this  limit,  which  recovers  an  earlier  result6,  is 
identical  to  the  condition  for  unstable  beam  envelope  oscillations,  derived 
for  such  beams  in  reference  7.  Note  that  we  consider  only  kq  >  0,  which 
indicates  right-handed  helicity  in  the  quadrupole  windings.  By  symmetry, 
(kq,B,)  is  equivalent  to  (-kq,-B_). 

*Work  supported  by  the  Defense  Advanced  Research  Projects  Agency. 
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Theoretical  results  suggest 
that  the  PoCE  can  operate  entirely 
in  the  stable  regime.  They  also 
suggest  that  as  y  is  increased, 
the  conditions  for  stability 
become  more  restrictive,  forcing 
operation  at  longer  pitch  lengths  ' '  f 

and  lower  quadrupole  gradients. ^  -T"  .5- 

This  may  not  present  a  problem  as  ‘=  . ; 

beam  focusing  requirements  are  CL  "[ 

generally  less  stringent  as  y  is  c- 

increased.  .  ■ 


NUMERICAL  RESULTS 

Simulations  have  been 
performed  for  parameters  typical 
of  the  SLIA  PoCE:  I ^  =  10  kA, 
normalized  RMS  emittance,  = 

.158  cm-rad,  r^  =  3  cm,  Bqk’  =  200 
G/cm,  It  =  0.1  and  Bz  =  5  kG.  Ve 
simulated  both  y  =  7  and  y  =  13, 
corresponding  to  the  beam  energies 
for  the  initial  and  final 
stellarator  sections  of  the  PoCE. 


Fig.  1  Stability  diagram  for 
beam  energy  y  =  7,  quadrupole 
gradient  Bqkq  =  200  G/cm  and 
waveguide  raaius  rg  =  3  cm. 


The  case  of  Bz  =  -5  kG  at  y  =  7, 
corresponding  to  the  smaller  (and 

less  useful)  stable  regime  in  Fig.  1,  was  also  simulated.  Stable  behavior 
was  observed  in  all  three  cases.  In  order  to  observe  transitions  between 
the  various  stability  regimes,  ve  set  kq  =  0.5  cm-1  and  y  =  7  and  considered 
Bz  =  1,  2,  4,  4.5  and  6  kG.  The  (kq,BzJ  parameter  space  locations  for  these 
runs  are  noted  in  Fig.  1.  The  results  are  listed  in  Table  1.  Shown  are 
theoretical  and  numerical  values  for  the  growth  rate  and  peak  frequency  for 


the  electromagnetic  TEjj  mode. 


Linear  Theory 


Simulation 


(kG)  T  (cm~M  ggeakZS  (cnri)  r  (cm-1)  “aeak£S_LHL-l 


peak 


1.0 

.0334 

0.71 

0.009 

0.75 

2.0 

.0438 

1.03 

0.019 

0.94 

4.0 

orbit-unstable 

- 

orbi t-unstabie 

- 

4.5 

.0266 

0.84 

0.020 

0.83 

6.0 

stable 

- 

stable 

- 

Table  1.  Theoretical  and  numerical  values  of  peak  growth  rates  and 
frequencies  for  varying  Bz. 
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Grovth  rates  were  measured 
by  analyzing  the  TE^  mode,  for 
which  the  B_  and  Er  components 
may  be  "projected  out"  from  the 
electromagnetic  spectrum  in  a 
straightforward  manner.  The  Bz 
component  of  the  TE^  node  is 
then  decomposed  via  FFT  into 
discrete  wavenumbers,  k.  The 
growth  rate  is  then  obtained  as  a 
function  of  k  for  those  cases 
where  an  unstable  mode  grew  above 
the  background  "noise". 

Similarly,  growth  rates  may  be 
obtained  as  a  function  of 
frequency,  o>.  Because  our 
simulation  takes  place  in  a 
coordinate  system,  (r,0, C=ct-z) , 
that  moves  with  the  beam,  the 
growth  rates  that  result  from 
this  analysis  of  the  data  are  T  = 
Im(oj/c  -  k) .  This  corresponds  to 
the  theoretical  result  only  in 
the  case  that  Im(k)  =  0  as  was 
assumed  in  obtaining  the 
theoretical  growth  rates  given  in 
Table  1. 


-4  - 


-10 


0.0  100  too  000  400  500 

Z  (cm) 

Fig.  2  The  Fourier  component 
of  at  u/c  =  0.83  cm~- 

plotteri  versus  z  for  the 
3,  =  4.5  kG  case. 


As  an  example,  the  numerically  measured  TS^i  mode  at  co/c  =  0.83  cm-1  is 
plotted  on  a  logarithmic  scale  versus  position,  z,  in  Fig.  2  for  the  B_  = 

4.5  kG  case.  From  this  plot,  the  mode  appears  to  saturate  at  z  =  350  cm. 
This  saturation  may  be  a  numerical  effect  due  to  the  fact  that  the  particle 
simulation  is  inherently  noisy.  A  carefully  scaled  parameter  study  would  be 
required  to  determine  actual  saturation  levels  for  the  instability. 


Several  aspects  of  the  simulation  results  deserve  comment: 

1)  To  minimize  computational  expense,  the  typical  simulation,  which 
takes  place  in  the  "beam  frame",  was  of  a  short  (<1  m)  beam  propagating  over 
a  distance  of  2-10  meters.  In  the  PoCE,  the  pulse  length  is  =  25  m  and  the 
interaction  length  (curved  section)  is  =  2.5  m. 

2)  Apparent  saturation  of  the  instability  was  observed  in  each  of  the 
three-wave  unstable  cases.  This  was  associated  with  a  beam  centroid 
displacement  in  the  1-2  mm  range,  an  emittance  increase  of  20-100%  and,  in 
most  cases,  the  onset  of  current  loss.  The  exact  saturation  mechanism, 
which  may  be  due  to  numerical  effects  or  to  selective  loss  of  beam 
particles,  was  not  determined. 

3)  The  spectrum  of  unstable  modes  was  typically  much  broader  than 
predicted  by  the  linear  theory.  This  was  especially  pronounced  in  the  Bz  = 

1  kG  and  B,  =  2  kG  cases,  where  the  growth  rates  were  significantly  lower 
than  predicted.  The  broad  spectra  are  not  surprising  in  that  the  beam 
radius,  r^  =  1-2  cm,  is  much  larger  than  is  assumed  in  the  theory. 
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4)  For  the  run  at  Bz  *  4  kG,  orbit  instability  vas  observed  as 
predicted  by  a  numerical  solution  of  the  dispersion  relation.  This  result 
disagrees  with  Fig.  1,  which  is  derived  in  the  zero  current  limit.  As  beam 
current  is  increased,  the  stability  boundaries  shift  in  such  a  way  that  the 
unstable  regions  become  larger.  At  1^  =  10  kA,  the  orbit  unstable  regime 
encompases  the  the  Bz  =  4  kG  point,  which  is  close  to  the  orbit  unstable 
boundary  at  zero  current.  In  this  case,  the  orbit  instability  caused  severe 
disruption  of  the  beam  envelope,  with  the  major  radius  of  the  beam  ellipse 
expanding  to  make  contact  with  the  wall  within  the  first  meter  of 
propagation. 

5)  Further  runs  at  PoCE  parameters,  y  <  13,  kq  =0.1,  Bz  =5  kG  and  100 
G/cm  <  Bqkq  <  200  G/cm,  were  stable  in  all  cases. 

CONCLUSIONS 

In  this  study,  we  have  used  electromagnetic  particle. simulation  to 
confirm  and  clarify  the  theoretical  results  of  references  2-4.  Ue  have  seen 
that  the  stability  boundaries  outlined  in  reference  4  are  good  predictors  of 
behavior  in  the  simulations.  More  importantly,  we  have  verified  that  the 
three-wave  stable  regime,  initially  discussed  in  reference  3,  lies  in  a 
region  of  parameter  space  suitable  for  the  PoCE.  This  regime  appears  at  Bz 
>  0  in  Fig.  1.  Note  that  with  our  sign  convention,  Bz  >  0  corresponds  to 
what  is  called  the  "reversed  field"  case  in  reference  3.  Two  minor 
departures  from  the  linear  theory  were  found.  Firstly,  unstable  cases 
generally  showed  lower  growth  rates  and  broader  spectra  than  expected. 
Secondly,  apparent  saturation,  associated  with  emittance  growth  and  the 
onset  of  beam  loss,  was  observed  in  each  three-wave  unstable  case.  Stable 
behavior  was  observed  in  both  stable  regimes.  Finally,  these  simulations 
confirm  theoretical  predictions^ » 4  that  the  SLIA  proof  of  concept  experiment 
can  operate  entirely  in  the  stable  regime. 
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THREE-DIMENSIONAL  SIMULATIONS  OF  BEAM  TRANSPORT  AND  EMITTANCE 


GROWTH  IN  THE  SPIRAL  LINE  INDUCTION  ACCELERATOR* 


J.  KRALL  and  G.  JOYCE 

Naval  Research  Laboratory 
Beam  Physics  Branch 
Plasma  Physics  Division 
Washington,  DC  20375-5000 

Magnetic  transport  of  high-current  electron  beams  is  a  problem  of 
increasing  importance  for  induction  accelerator  physics.  An  example  is  the 
spiral  line  induction  accelerator  (SLIA)l,  which  uses  solenoidal  focusing  in 
the  straight  sections  and  both  solenoidal  and  stellarator  (helical 
quadrupole)  focusing  in  the  curved  sections.  Potential  difficulties  with 
this  scheme  arise  from  the  need  to  transport  the  beam  through  transitions 
from  solenoidal  to  stellarator  focusing  on  each  turn.  The  need  here  is  to 
minimize  or  eliminate  emittance  growth  through  each  transition. 

The  present  numerical  study  is  the  product  of  ELBA^ ,  a  three-dimensional 
particle  code  recently  developed  at  NRL.  ELBA  simulates  a  beam  propagating 
within  a  cylindrical  metallic  pipe.  The  full  set  of  Maxwell's  equations  and 
the  full  relativistic  motion  of  the  beam  particles  are  included.  The  beam  is 
matched  to  the  field  configuration  to  minimize  initial  oscillations .  This 
matched  configuration  is  calculated  by  STELMAT^,  which  was  adapted  for  use  in 
ELBA. 


BEAM  TRANSPORT  IN  STELLARATOR  FIELDS 


Previous  studies  of  beam  transport  in  stellarator-focused  accelerators 
have  been  primarily  concerned  with  stellarator-betatron  configurations,  where 
transitions  from  stellarator  to  solenoidal  focusing  are  not  required  as  in 
the  SLIA.  This  matching  problem  is  ideally  solved  through  the  use  of 
focusing  elements  which  somehow  convert  a  matched  beam  in  solenoidal  fields 
(circular  cross  section)  into  a  matched  beam  in  the  stellarator  fields 
(elliptical  cross  section).  Recent  studies^  have  suggested  strategies 
through  which  this  may  be  accomplished,  but  without  directly  addressing  the 
issue  of  emittance  growth  through  these  transitions.  Intuitively,  we  expect 
configurations  for  which  the  matched  beam  ellipse  has  a  high  eccentricity  to 
cause  the  greatest  emittance  growth. 


In  this  preliminary  study,  we  use  the  ELBA  simulation  code  to  transport 
a  short  slice  of  beam  through  stellarator  fields  typical  of  those  discussed 
for  the  SLIA  proof  of  concept  experiment  (PoCE).*’4  Two  significant 
approximations  are  used.  Firstly  we  ignore  curvature  effects  by  performing 
the  simulations  in  straight  cylindrical  geometry.  Secondly  we  idealize  the 
external  fields,  ignoring  fringing  fields  and  other  realistic  corrections. 

In  each  case  quadrupole  fields  are  ramped  smoothly  with  an  envelope  given  bv 


*Uork  supported  by  the  Defense  Advanced  Research  Projects  Agency. 
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(1) 


«»>  -  H1  * 

ramp 

where  the  entrance  to  the  stellarator  section  is  located  at  z  =  Z]  and  zramp 
is  a  scale  length  over  which  this  function  increases  from  0.12  to  0.88.  The 
quadrupole  fields  are  "ramped  down"  in  a  similar  fashion  at  the  exit.  In  all 
cases  we  chose  zramp  =  6  cm,  and  fixed  the  length  of  the  stellarator  section 
at  aR  -  "’51.33  cm,  where  R  =  80  cm  is  the  radius  of  curvature  of  the  PoCE 
bends . 


In  this  study,  we  will  use  an  emittance  definition 
coupled  x-y  motion^: 


for  beams  with 


e 

xy ,  rms 


where  the  matrix  I  is  given  by 


A/  H 


I.  .  =  <u .  u  •>  -  <u . Xu  .> 
ij  i  J  i  I 


(2) 


(3) 


Here,  "<>"  signifies  an  average  over  the  beam  particles  and  i,j  =  1,2, 3, 4 
with  uj  =  x,  U2  =  dx/dz,  U3  =  y  and  U4  =  dy/dz.  The  normalized,  RMS 
emittance  is  £n>rms  =  3r£Xy  ,  rms  • 

NUMERICAL  RESULTS 


We  first  consider  a  high  eccentricity  case.  The  parameters  are  beam 
current  1^  =  10  kA,  relativistic  factor  y  =  7,  en  rms  =  0.158  cm-rad,  axial 
field  Bz  =  5  kG,  quadrupole  gradient  V3q  =  198  G/cm  and  quadrupole  pitch 
length  Ap  =  125.7  cm.  In  this  case  the  matched  beam  parameters  in  the 
stellarator  section  are  Rmajor  =  1.796  cm  and  Rminor  =  0.514  cm  so  that  the 
eccentricity,  given  by 


a 


1 


R~.  , 

minor | 

R2  .  ' 

major 


1/2 


(4) 


is  a  =  0.958.  The  major  and  minor  radii  for  this  case  are  plotted  versus  z 
(propagation  distance)  in  Fig.  1.  The  normalized  emittance,  plotted  versus  z 
in  Fig.  2,  shows  a  47%  increase. 


The  equilibrium  beam  eccentricity  for  these  parameters  can  be  lowered  by 
decreasing  either  the  pitch  length  or  the  quadrupole  gradient.  The  emittance 
increase  is  plotted  versus  pitch  length  in  Fig.  3.  The  eccentricity,  a,  is 
also  given.  Plots  versus  quadrupole  gradient  are  shown  in  Fig.  4.  In  both 
cases,  the  emittance  variation  falls  below  measurable  levels  (<3%)  at 
a  =  0.65. 
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Low  eccentricity  stellarator  configurations  have  also  been  studied  for  the 
PoCE.  A  set  of  such  parameters  for  each  of  the  three  PoCE  bends  are  given  in 
Table  1.  Simulation  results  showed  that  a  y  =  7,  1^  =  10  kA,  cn  rms  =  0.158 
cm-rad  beam  was  transported  through  the  first  bend  with  no  measurable 
emittance  growth.  In  fact,  as  the  quadrupole  gradient  was  increased  from 
115.7  to  350  G/cm  in  this  case  (a  increasing  from  0.507  to  0.835),  the 
emittance  growth  remained  negligible.  Note  that  the  quadrupole  field  period, 
Xp/2  =  15  cm,  is  not  significantly  larger  than  the  ramp-up  length  for  the 
quadrupole  fields  (zramp  =  6  cm)  in  this  case.  This  suggests  that  higher 
eccentricity  stellarators  are  tolerable  when  0(Xp/2)  =  0(zramp). 


Bend 

Y 

Bz  (kG) 

VB^  (G/cm) 

Xp  (cm) 

^maj 

^min  (cm) 

a 

1 

7 

4.666 

115.7 

31.42 

.776 

.669 

.507 

2 

10 

4.228 

103.5 

50.27 

.851 

.677 

.606 

3 

13 

4.446 

107.6 

62.83 

.854 

.645 

.655 

Table  1.  Low  eccentricity  parameters  for  each  of  the  three  PoCE  bends. ^ 


Finally  we  performed  a  simulation  of  beam  transport  through  three  PoCE 
bends,  using  parameters  given  in  Table  1,  with  y  =  7,  1^  =  10  kA  and 
cn  rms  =  0-158  cm-rad.  The  beam,  which  was  accelerated  between  each 
stellarator  section,  showed  no  measurable  emittance  growth. 

CONCLUSIONS 

This  preliminary  study  suggests  that  for  poorly  chosen  parameters, 
transport  through  the  stellarator  focusing  sections  can  cause  emittance 
growth  as  high  as  47%.  Such  parameters  feature  high  eccentricity  beam 
equilibria,  a  >  0.65,  and  long  pitch  lengths,  Xp/2  >>  zramp>  where  zramp  is 
the  ramp-up  length  for  the  quadrupole  fields.  Vhen  the  quadrupole  fiexd 
period,  Xp/2,  is  of  the  order  of  the  ramp-up  length  for  the  quadrupole 
fields,  high  eccentricity  (a  >  0.8)  stellarator  sections  are  tolerable. 
Successive  transitions  through  three  low-eccentricity  stellarator  sections 
showed  negligible  emittance  growth.  It  is  important  to  keep  in  mind, 
however,  that  the  external  fields  and  geometry  used  in  this  study  were  highly 
idealized . 
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NRL  STUDIES  OF  RLA  ACCELERATOR  PHYSICS* 

Glenn  Joyce,  Richard  Hubbard,  and  Richard  Fernsler 

Beam  Physics  Branch,  Plasma  Physics  Division 
Naval  Research  Laboratory,  Washington,  D.C. 

INTRODUCTION 

One  proposed  design  for  the  Recirculating  Linear  Accelerator  is  a  race 

track  configuration.  An  IFR  channel  is  initially  set  up  in  the  entire  device 

using  a  Godfrey  gun.  The  channel  must  remain  intact  throughout  the  entire 

experiment  which  lasts  about  400  ns  and  consists  of  four  transits  of  the  beam 

around  the  racetrack.  The  behavior  of  the  ion  channel  over  this  time  period 

is  of  some  concern  since  the  channel  will  help  guide  the  beam  and  provide 

space  charge  cancellation.  A  simple  calculation  of  the  time  for  the  radius  to 

2 

double  due  to  its  own  self  forces  is  t  -  (a  /vf  m./m  )l/2  -  44  ns.  This 

i  e 

estimate  is  pessimistic  since  it  assumes  that  the  channel  will  maintain  its 
shape  during  the  acceleration  time.  In  fact,  the  outer  ions  will  travel  at 
this  rate,  but  the  inner  ions  will  hardly  move  at  all  so  it  is  unlikely  that 
the  channel  will  maintain  its  shape.  The  behavior  of  the  interior  of  the 
channel  is  probably  more  important  than  that  of  the  channel  edge. 

CHANNEL  NEUTRALIZATION 

If  the  channel  becomes  neutralized  after  the  passage  of  the  beam,  the  time 
scale  will  become  much  longer,  but  there  is  a  question  of  where  the 
neutralizing  electrons  will  come  from.  The  best  case  scenario  is  that  all  the 
channel  electrons  will  go  to  the  chamber  walls  and  remain  there  as  the  beam 
passes  by.  After  the  beam  passes,  the  channel  ions  will  pull  the  electrons 
back  toward  the  channel  center.  However,  the  potential  difference  for  a  10 
ka.  beam  is  -  230  Kev  and  the  kinetic  energy  of  these  electrons  will  be  at  its 
maximum  when  they  are  in  the  center  of  the  channel  so  there  will  be  little 
tendency  for  them  to  stop  on  the  channel. 

In  order  to  investigate  the  behavior  of  the  channel  for  long  times,  we 
have  done  a  number  of  numerical  simulations  using  the  FRIEZR  code,  an 
axisymmetric,  relativistic  particle  simulation  code.  We  have  performed  a  set 
of  scaled  "experiments"  using  the  ion  mass  of  Hydrogen  instead  of  that  of  a 
heavier  mass  such  as  Xenon.  Since  the  time  scales  vary  as  the  square  root  of 
the  mass  ratio,  we  have  scaled  back  all  the  times  by  a  factor  of  ten.  As  a 
result,  the  beam  length,  its  erosion,  transport,  and  other  properties  will  be 
incorrectly  treated,  but  the  behavior  of  the  channel  ions  should  give  us  an 
indication  of  what  will  happen  in  the  actual  experimental  device. 
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Ve  have  chosen  the  following  set  of  numerical  parameters: 


I  = 

8  kA 

Y  =  10 

ab  - 

1  m 

a  ,  =  1  cm 
ch 

f  = 

.3 

m./m  =  2000 

.  /  in 

1  e 


1,  =  2.5  ns 

D 

which  corresponds  to  an  "equivalent"  set  of  experimental  parameters 

m./m  =  ( 100)*2000  1.  =  25  ns 

l  e  b 

pulse  separation  =  75  ns. 


The  simulation  results  are  summarized  in  Figs.  1  and  2.  Over  the  time  of 
two  transits  of  the  beam,  the  wings  of  the  channel  have  spread  out 
significantly,  leaving  behind  a  weak  and  narrow  central  core  which  will  slowly 
decay  away.  The  electrons  after  two  transits  have  spread  out  almost  uniformly 
across  the  chamber  with  little  indication  of  a  higher  density  at  the  channel. 
For  our  parameters,  this  represents  about  a  1CT  neutralization  of  the  channel. 
It  appears  that  without  some  strategy  for  replenishing  the  ion  channel,  the 
ion  guiding  effects  may  become  too  weak  to  be  effective  in  providing  space 
charge  neutralization  and  guiding  for  the  near  term  exper ime.n ts . 


o.;t 


iuii  oeiijirr 


o.i 


0.00 


0.06 


0.01  i- 

I 


0  r 
-o 


1 1 
'  I 


0,1  r  EltClIiUtl  UCNSITT 


I  C 


Ion  density  vs.  radius 


r.'  ~ 

-  *  a  ■ 


Electron  density  vs.  radius 


22 


In  both  figures,  the  dashed  line  is  the  density  at  t  =  0,  and  the  solid  line 
is  the  density  at  t  =  200  ns. 

CHANNEL  STRENGTHENING 

It  is  possible  that  the  beam,  itself,  may  be  used  to  strengthen  the  channel 
during  each  of  its  passes  around  the  racetrack..  One  strategy  is  to  provide  a 
neutral  density  high  enough  that  the  beam  will  ionize  it  to  some  extent 
generating  an  enhanced  channel  at  each  pass.  Beam  ionization  would  require 
neutral  densities  much  higher  than  have  been  planned.  If,  for  example,  we 
allow  the  beam  to  ionize  up  to  f  =  .2,  the  background  density  must  be  -  8 
microns.  Another  possibility,  suggested  by  Sandia,  is  to  shorten  the  length 
of  the  track  so  that  the  time  between  passes  is  about  25  ns  so  that  channel 
will  not  have  expanded  significantly  before  the  beam  returns.  In  this  case, 
the  attraction  of  the  negatively  charged  beam  will  reverse  the  expansion  of 
the  channel.  Ve  have  done  simulations  for  both  cases,  and  in  both  cases,  the 
channel  ion  density  will  be  increased  to  provide  space  charge  neutralization. 
However,  in  neither. case  will  the  channel  guide  the  beam  to  the  center  of  the 
transport  vessel  if  the  channel  has  drifted  off  center  between  passes. 

Figures  3  and  4  show  surface  plots  of  the  ion  density  as  a  function  of 
time.  In  Fig.  3,  there  is  no  significant  beam  ionization,  while  in  Fig.  4, 
the  neutral  gas  density  is  high  enough  for  the  beam  to  generate  a  channel  of 
strength  f  =  .2. 


I OH • CHANNEL  DENSITY 


ION  CHANNEL  DENSITY  WITH  BEAN  GENERATED  IONS 
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CHANNEL  DRIFT 

A  major  concern  remains  about  the  drifting  of  the  channel  centroid  during 
the  experiment.  Although  ion  hose  may  not  be  a  problem  for  beams  vhose  length 
is  -  25  ns,  the  beam  and  ion  channel  do  interact.  Simulations  by  John  Vagner 
and  Leon  Feinstein  have  indicated  that  there  is  some  channel  drifting  after 
the  passage  of  the  beam.  There  is  no  restoring  force  to  center  the  channel  in 
the  beam  line.  In  fact,  if  there  is  not  a  significant  number  of  electrons, 
the  electrostatic  image  forces  due  to  the  bare  ion  channel  tend  to  accelerate 
it  tovard  the  vail.  We  can  get  an  estimate  of  thetime  for  the  channel  to  hit 
the  wall.  The  result  is 

cx  -  (I./2I  m./m  ) 1/2  b  ln{(yo/b  +  v(I  /2I  m./m  )  1  / 2 ) - 1 ) 

A  C  1  6  AC-k.6 

where  I  =  f  I, 
c  b 

I,  =  17  kA 
b  =  vail  radius 
yo  =  initial  offset 
y  =  initial  velocity 

An  estimate  of  the  time  for  our  parameters  is  x  -  250  ns. 

The  consideration  of  channel  motion  vill  require  mere  extensive  investigation. 
CONCLUSIONS 

Our  analysis  of  ion  channel  motion  in  the  Recirculating  Linear  Accelerator 
has  shown  that  electrons  from  the  chamber  vails  vill  not  be  effective  in 
reducing  the  spread  of  the  ion  channel  due  to  the  electrostatic  repulsion  of 
the  ions  although  there  may  be  methods  of  replenishing  the  ions  or  reducing 
the  spreading  by  using  tr.e  beam  itself.  The  problem  of  ion  channel  drift  may 
be  significant  although  more  study  vill  be  necessary. 

*Work  supported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA  Order 
No.  A395,  Amendment  80,  monitored  by  Naval  Surface  Warfare  Center. 
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EROSION  OF  ELECTRON  BEAMS  IN  THE  RLA  ACCELERATOR* 

R.  F.  Hubbard,  G.  Joyce  and  S.  P.  Slinker 
PLasma  Physics  Division,  Naval  Research  Laboratory,  Washington,  DC  20375 

I.  INTRODUCTION 

The  Recirculating  Linear  Accelerator  (RLA)  currently  being  developed  at 
Sandia  National  Laboratories  will  utilize  ion-focused  regime  (IFR)  transport 
in  both  straight  and  curved  sections  of  the  accelerator.  In  the  curved 
sections,  a  vertical  field  approximately  matched  to  the  nominal  beam  energy 
and  major  radius  of  curvature  of  the  bend  will  augment  the  IFR  channel.  In 
this  paper,  we  describe  ELBA  simulation  studies  of  beam  transport  and  erosion 
in  the  RLA.  The  purposes  of  these  studies  are  to  determine  beam  head  erosion 
rates  for  the  RLA  in  both  straight  and  curved  sections,  to  estimate  beam 
drifts  in  bends  arising  from  such  sources  as  local  vertical  field  mismatch  and 
hoop  forces,  to  estimate  emittance  growth  during  transport  and  to  study  beam 
degradation  in  the  broad,  weak  channels  expected  after  the  beam  recirculates. 

II.  METHODOLOGY 

ELBA  simulation  code  description:  ELBA  is  a  3-D  version  of  the  FRIEZR 
axisymmetric  particle  simulation  code.  It  is  fully  electromagnetic  and  solves 
the  field  equations  in  r,  0,  C  =  ct  -  z  and  t.  This  doppler-shif ted 
coordinate  system  is  particularly  well  suited  to  treating  IFR  transport  since 
it  moves  at  the  speed  of  light,  making  the  coordinate  C  an  approximate 
constant  of  the  motion  for  beam  electrons.  Although  z  is  not  a  dynamical 
variable,  diagnostics  are  usually  transformed  to  the  frozen  (z,  C)  coordinate 
system  to  facilitate  comparison  with  experiments.  The  field  solver  in  curved 
sections  assumes  a  small  aspect  radio  r/R  where  r  <  a  is  the  distance  from 

the  axis,  a  is  the  pipe  radius  and  R  is  Hie  mainr  radius  *■>(  the  bend.  ELBA 

w  1 
employs  an  elegant  mapping  algorithm  de-eloped  by  Friedman  for  treating  beam 

dynamics  in  the  bends.  The  implementation  of  this  method  in  ELBA  is  described 

in  more  detail  in  Ref.  2.  Beam  and  plasma  electron  dynamics  are  treated  fully 

relativistically ,  and  ion  motion  is  neglected  since  we  are  concentrating  on 

3 

beam  head  effects.  (Ion  channel  motion  studies  are  reported  elsewhere).  The 
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ratio  of  the  bending  field  to-  its  matched  value  =  By^mcVeR  can  be 
specified  arbitrarily. 

Erosion  Rate  Diagnostics:  In  the  RLA,  the  vail  radius  is  not  much  larger 

than  the  beam  radius  a^,  so  the  full  vidth  half  maximum  (FVHM)  of  the  beam 

current  I^(C,  2)  is  a  convenient  diagnostic  for  calculating  the  simulation 

erosion  rate  0  .  If  C  (2)  is  defined  by  I ( C  ,2)  =  I  / 2 ,  then 
e  p  p  max 

0e  =  “  ^p(zl>)/(z2  ~  Zl^  _  av/R>  ^ast  term  adiusts  f°r  slippage 

of  the  coordinate  system  in  bends. ^  Results  can  be  compared  vith  the 

theoretical  inductive  erosion  rate:  0  (th)  =  fLv/y-  Here  f  is  the  ion-to- 

e 

beam  density  ratio,  v  =  I,/ (17  kA)  and  the  inductance  L  =  21n(a  /a. ). 

b  v  b 


III.  SIMULATION  RESULTS 


Beam  and  channel  parameters:  A  series  of  simulations  were  run  vith 

Yn  =  10,  a,  =  1  cm,  a  =  4.5  cm,  R  =  84  cm,  I  =  10  kA  and  channel  radius 
U  b  v  max 

ach  "  1  cm.  A  128  cm  long  beam  vas  transported  500  cm  in  straight  or  curved 
sections.  The  500  cm  in  a  curve  corresponds  to  340°  although  a  single  bend 
in  the  experiment  is  only  180°  Results  from  the  various  simulations  are 
summarized  in  the  table  belov. 


Transport 


Case 

f 

B  /B 
y  m 

0e 

Efficiency 

Cr 

Comments 

FI 

0.5 

N.  A. 

0.079 

>0.99 

5 

Straight  section,  0  (th)=0.079 

F2 

0.5 

1.0 

0.119 

>0.99 

5 

Curved;  no  sector  magnet  corr. 

F3 

0.5 

1.0 

0.110 

>0.99 

5 

Curved;  sector  magnet  on 

F4 

0.25 

1.0 

0.063 

>0.99 

5 

Curved;  lover  f 

F5 

0.5 

1.1 

0.121 

0.99 

5 

Curved;  stronger  vertical  field 

F6 

0.25 

1.1 

0.063 

0.99 

t; 

Curved;  lover  f  and  higher  B 

F7 

0.5 

0.9 

0.106 

0.99 

5 

Curved:  veaker  vertical  field 

SI 

0.25 

1.0 

0.041 

0.96 

60 

Curved:  std.  moderate  rise  case 

S2 

0.25 

N.  A . 

0.030 

0.95 

60 

Straight;  0  (th)  =  0.040 

S3 

0.25 

0.9 

0.042 

0.89 

60 

Curved;  veafer  vertical  field 

S4 

0.25 

1.1 

0.059 

0.96 

60 

Curved:  stronger  vertical  field 

VI a 

0.10 

N.  A. 

0.010 

0.84 

60 

Straight;  0  (th)  =  0.008 

V2a 

0.10 

1.0 

0.013 

0.84 

60 

Curved;  std.  veak  channel  case 

V3a 

0.10 

1.0 

0.014 

0.69 

60 

Curved;  assumes  407.'  inc.  in  z 

■  3. 

0.10 

1.0 

0.030 

0.48 

60 

Curved;  Yq  =  10 

Curved;  veaker  vertical  field 

V6a 

0.10 

0.9 

0.022 

0.29 

60 

V  6a 

0.10 

1.1 

0.007 

0.36 

60 

Curved;  stronger  vertical  field 

aVeak 

channel  "recirculated"  beam  case:  a  . 

=  3.5  cm,  Yq  =  21:  normalized 
aas  in  previous  cases). 

emittance  z 

matched  for 
n 

f=0.25  and 

O 

M 

►— * 
O 

r 
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Erosion  results  for  fast  rise  cases:  In  order  to  facilitate  comparison 

with  analytical  results,  a  series  of  simulations  with  a  very  short  (5  cm) 

current  rise  length  were  carried  out.  For  a  channel-to-beam  density  ratio 

f  =  0.5  (Case  FI),  the  erosion  rate  in  a  straight  section  simulation  agreed 

with  the  analytical  estimate  Bg  =  0.079.  The  rates  in  the  curved  section 

f  =  0.5  simulations  (Cases  F2,F3  F5  and  F7)  were  typically  40X  higher  than  in 

the  straight  case.  As  expected,  the  erosion  rate  dropped  by  almost  a  factor 

of  two  when  f  was  reduced  to  0.25  (Cases  F4  and  F6).  Some  cases  included  a 

10%  mismatch  between  the  vertical  field  B  and  the  matched  value  B  . 

y  m 

Erosion  and  beam  loss  in  moderate  rise-rate  cases:  These  cases  (S1-S4) 

are  similar  to  those  above  except  that  Cr  =  60  cm.  Erosion  rates  are  somewhat 

lower  than  those  produced  in  the  previous  cases.  Figure  1  plots  the 

displacements  X(C)  and  Y(0  for  a  matched  bending  field  simulation  (Case  SI). 

Ohmic  energy  loss  cases  the  beam  head  to  drift  inward  (X  <  0)  while 

displacements  in  Y  remain  small.  Ve  believe  that  the  small  positive  x- 

deflection  in  the  beam  body  seen  in  Fig.  1  is  due  to  the  hoop  force.  Choosing 

B  /B  =  0.9  instead  of  1.0  (Case  S3)  results  in  a  smaller  inward  drift  in  the 
y  m 

beam  head  and  a  larger  outward  drift  in  the  body. 


Fig  1.  Beam  displace¬ 
ment  X(C)  (solid)  and 
Y(C)  (dashed)  at 
z=300cm:  Case  SI 


Fig  2.  Pinched  beam 
and  escaping  plasma 
electrons  at  z=318cm 
and  C=64cm:  Case  SI 


Fig  3.  Beam  displacement 

X(C)  and  Y(z)  at  z=300cm 

for  Case  V2  (Veak  ion 

channel;  B  =B  ) . 

y  m 


Figure  2  plots  the  position  of  simulation  particles  at  z 


313  cm  and  C  = 


64  cm.  Beam  electrons  are  mostly  confined  to  the  area  near  the  center,  while 


the  vertical  "column"  of  particles  consists  of  plasma  electrons  escaping  along 


the  vertical  field  lines. 


Erosion  and  beam  loss  in  a  weak  ion  channel:  Ion  channel  expansion 
between  recirculations  is  a  major  concern  for  the  racetrack  design.  Cases  VI- 


>, 
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V6  mimic  this  effect  by  using  a  higher  energy  beam  and  a  broader,  weaker  ion 
channel.  As  expected,  these  cases  have  substantially  lower  erosion  rates  but 
more  particle  loss  to  the  wall.  Particle  loss  is  primarily  due  to  evaporation 
and  drifting;  it  is  particularly  severe  in  Cases  V4  (low  beam  energy)  and  V 5 
(low  vertical  field).  In  addition,  the  ion  channel  is  much  less  effective  in 
centering  the  beam,  as  can  be  seen  in  the  matched  bending  field  case  V2  in 
which  the  beam  centroid  is  displaced  by  more  than  1  cm  (Fig.  3).  This  outward 
drift  is  due  both  to  the  hoop  force  and  to  a  small  increase  on  beam  energy 
arising  from  inductive  fields  produced  by  rapid  beam  expansion  just  after  the 
beam  is  injected. 

There  is  very  little  growth  in  rms  emittance  even  though  the  beam  is  not 
well  matched  to  the  channel  strength.  Particles  which  would  normally 
contribute  to  emittance  growth  are  generally  lost  to  the  wall  instead.  The 

4 

results  are  consistent  with  end-to-end  BUCKSHOT  simulations  by  Vagner. 

IV.  CONCLUSIONS  AND  REFERENCES 

Simulations  of  the  RLA  using  EL3A  generally  produce  erosion  rates  in 
the  straight  sections  which  agree  with  the  standard  inductive  erosion  formula. 
In  curved  sections  is  generally  30-50"  higher  than  in  straight  sections.  A 
10"  mismatch  in  cases  a  small  shift  in  the  beam  centroid  and  has  a  modest 
effect  on  g  .  Weak  channel  ("recirculated")  beam  simulations  show  little 
erosion  or  emittance  growth,  but  often  exhibit  large  drifts  from  the  axis  and 
substantial  current  loss. 

1.  A.  Friedman,  D.  Grote  and  I.  Haber,  Bull.  Am.  Phys.  Soc.  34,  2064  (1989). 

2.  G.  Joyce,  et  al.,  "ELBA  -  a  Three-Dimensional  Particle  Simulation  Code," 
these  proceedings. 

3.  G.  Joyce,  et  al.,  "IFR  Propagation  of  Beams  as  Applied  to  RLA,"  these 
proceedings . 

4.  J.  Vagner,  "Simulations  for  RLA  Experiment  ana  Prototype,",  these 
proceedings . 

Vork  supported  by  the  Defense  Advanced  Research  Projects  Agency.  ARPA  Order 
No.  4295,  Amendment  No.  80,  monitored  by  the  Naval  Surface  Varfare  Center. 
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ATOMIC  AND  PLASMA  SURFACE  INTERACTIONS 
IN  ACCELERATORS  AND  IFR  CELLS* 


by 

A.  W.  ALI 


Naval  Research  Laboratory,  Washington,  DC  20375 


1.  INTRODUCTION 

DARPA  is  proceeding  with  its  support  and  funding  of  new 
concepts  in  compact  accelerators.  Clearly  the  dominant  issues  for 
accelerators  fall  within  the  domain  of  pure  plasma  physics  where 
electrons  are  moving  in  a  vacuum,  confined,  guided,  by  magnetic 
fields  and  accelerated  by  electric  fields.  However,  there  is  no 
such  thing  as  a  vacuum.  Beam  electrons  are  accelerated  in 
chambers  generally  evacuated  to  some  low  pressures  in  the  range 
of  10-5  -10-9  Torr.  Under  these  conditions,  atmospheric  species 
like  N2,  02  and  contaminants  as  H20,  C  and  other  gases  are 
present,  albeit  in  low  densities  in  the  range  of  107  -  1011  cm-3. 
These  gaseous  species  will  be  ionized  and  dissociated  by  the 
circulating  beam  electrons  resulting  in  a  plasma  (electrons  and 
positive  ions)  which  may  impact  the  beam  electrons,  their 
trajectories,  and  their  accelerations.  However,  there  is  another 
and  very  interesting  phenomena  which  involves  the  inner  surfaces 
of  the  accelerating  chambers  through  which  the  beam  electrons 

*This  work  was  supported  by  Defense  Advanced  Research  Project 
Agency  under  ARPA  order  No.  4395,  Amendment  No.  80  and  monitored 
by  the  Naval  Surface  Weapon  Center. 
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circulate.  This  phoenomena  is  equally  relevant  to  IRF  cells  which 
are  often  used  in  beam  conditioning  in  electron  beam  propagation 
experiments.  These  surfaces  act  as  a  source  of  electrons,  ions, 
and  neutrals  depending  on  various  atomic  processes  which  occur 
when  particles  and  radiation  fall  on  these  surfaces.  These 
surfaces  generally  acquire  a  monoatomic  layer  of  atoms  and 
molecules  which  are  adsorbed  into  the  surface.  The  adsorption 
occurs  in  time  as  a  result  of  the  gas  flow  in  and  out  of  the 
acceleration  chambers  and  IFR  cells  (filling  and  evacuation)  and 
by  ions  arriving  at  the  surface.  These  ions  acquire  an  electron 
and  are  neutralized.  Hence,  in  a  state  of  equilibrium  one 
envisions  a  neutral  species  density  distribution  being  maximum  at 
the  surface,  decaying  exponentially  toward  the  center  of  the 
chamber.  This  gaseous  distribution  is  what  one  calls  the 
outgassing  of  the  surfaces.  Depending  on  the  magnitude  of  the 
outgassing  from  the  surface,  considerable  ionization  can  result 
by  the  interaction  of  the  beam  electrons  with  these  species  which 
are  defusing  into  the  chamber. 

Additionally,  ions  accelerated  by  the  beam  will  charge 
exchange  with  the  neutrals  generating  fast  neutrals  which  strike 
the  chamber  surface  releasing  additional  atomic  species  into  the 
chamber.  Furthermore,  the  accelerated  electrons  will  emit 
radiation  which  upon  interaction  with  the  surfaces  generate 
additional  electrons  and  ions.  So  will  the  accelerated  plasma 
electrons.  Accordingly,  there  is  a  host  of  atomic  processes  which 
occur  in  any  accelerator  chamber,  IFR,  and  conditioning  cells 
which  must  be  considered  and  its  impact  on  beam  breakup  and 
breakdown  in  IFR  cells  quantified. 

2.  BASIC  ISSUES 

The  basic  physical  processes  which  lead  to  the  enhancement  of 
the  ionization  level  in  accelerators  and  IFR  cells,  change  the 
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channel  conductivity,  and  lead  to  the  beam  breakup,  require 
detailed  analysis.  The  physical  processes  are: 

1.  The  beam  electrons  colliding  with  the  ambient  gas  and 

the  outgassed  species  resulting  in  the  ionization  of 
the  neutral  species,  thus  generating  a  plasma. 

2.  Accelerated  plasma  electrons  and  ions  impacting  the 
chamber  surface,  resulting  in  the  ejection  of  atoms, 
negative  ions,  and  molecules  from  species  adsorbed  to 
the  surface.  Such  collisions  contribute  to 
adsorption,  desorption  or  sputtering,  reflection  and 
surface  heating.  These  processes,  in  general, 
increase  the  gas  density  near  the  surface  and  cause 
increased  ionization  by  the  collisions  of  charged 
particles  with  these  species.  The  density  of  the 
adsorbed  atoms  or  molecules  generally  is  in  the  order 
of  1015  cm-3,  based  on  a  monoatomic  layer  adsorbed  on 
the  surface.  The  gas  density  distribution  starts  from 
this  value  and  decays  exponentially  to  zero  (in 
perfect  vacuum)  at  the  center  of  the  accelerator  and 
IFR  chambers 

3.  Fast  neutral  species  impacting  the  electrode  surface. 
These  fast  neutrals  generally  arise  as  a  result  of 
charge  exchange  of  accelerated  ions  with  atoms  or 
molecules,  whereupon  the  fast  ion  becomes  a  fast 
neutral . 
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To  determine  the  effects  of  these  interactions  on  the 


beam  stability,  one  may  use  analytic  and  ab-initio  methods. 
The  ab-inito  approach,  in  general,  is  complicated  and 
requires  the  time  and  space  dependent  simulations  of 
electron  and  ion  orbits  in  the  accelerating  chamber  and 
their  interactions  with  gas  species  and  surfaces.  This  must 
be  coupled  to  basic  atomic  processes  in  order  to  describe 
the  ionization  and  the  breakdown  phenomenon.  The  simulation 
should  be  tailored  to  all  appropriate  accelerating  chambers 
and  IFR  cells.  In  addition  to  detailed  simulation  method, 
simple  analytic  methods1  can  be  utilized  to  obtain 
reasonable  estimates  for  ionization  rates  and  breakdown 
times  for  a  wide  range  of  electric  fields,  outgassed  species 
densities ,  chamber  diameter,  beam  circulation  time  and 
chamber  radius,  etc. 

1.  A.  W.  Ali,  Laser  and  Particle  Beams,  6,  105  (1988) 
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Advanced  Scientific  Concepts,  Inc. 

2441  Foothill  Lane 
Santa  Barbara,  CA  93105 
Attn:  Dr.  Roger  Stettner 

Advanced  Technologies  Research 
14900  Sveitzer  Lane 
Laurel,  MD  20707 
Attn:  Mr. 'Daniel  Weldman 

The  Aerospace  Corporation 
Mail  Stop  M2-269 
P.  0.  Box  92957 
Los  Angeles,  CA  90009 
Attn:  Dr.  David  L.  McKenzie 
Dr.  Carl  J.  Rice 

AFATL/DLJW 

Elgin  Force  Base,  FL  32542 
Attn:  MAJ  Louis  W.  Seller,  Jr. 

Air  Force  Office  of  Scientific  Research 
Physical  and  Geophysical  Sciences 
Bolling  Air  Force  Base 
Washington,  DC  20332 
Attn:  Major  Bruce  Smith 

Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117-6008 
Attn:  Dr.  William  L.  Baker  (AFWL/NTYP) 

Dr.  Brendan  B.  Godfrey 
Dr.  Inara  Kuck 

Applied  Physics  Laboratory 

The  Johns  Hopkins  University 

Asst,  to  Dir.  for  Tech.  Assessment 

Johns  Hopkins  Road 

Laurel,  MD  20707 

Attn:  Dr.  Samuel  Koslov 

Armed  Forces  Radiobiology 
Research  Institute 
Chief,  MRAD 
NMC-NCR 

8ethesda,  MD  20814-5145 
Attn:  LCDR  J.  P.  Jacobus 

U.  S.  Army  Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground,  Maryland  21005 
Attn:  Dr.  Donald  Eccleshall  (DRXBP-BM) 
Dr.  Anand  Prakash 
Dr.  Clinton  Hollandsvorth 


Avco  Everett  Research  Laboratory 
2385  Revere  Beach  Pkvy 
Everett,  Massachusetts  02149 
Attn:  Dr.  R.  Patrick 

Dr.  Dennis  Reilly 

Ballena  Systems  Corporation 
P.  0.  Box  752 
Alameda,  CA  94501 
Attn:  Dr.  Adrian  C.  Smith 

Dr.  William  E.  Wright 

Ballistic  Missile  Def.  Ad.  Tech.  Ctr. 
P.0.  Box  1500 

Huntsville,  Alabama  35807 
Attn:  Dr.  M.  Hawie  ( BMDSATC- 1 ) 

Dr.  M.  J.  Lavan  (BMDATC-E) 

Mr.  Dan  Whitener 

The  Boeing  Aerospace  Company 

MS-2E30 

Box  3999 

Seattle,  WA  98124 

Attn:  Dr.  Robert  C.  Milnor 

Booz,  Allen,  and  Hamilton 
Crystal  Square  2,  Suite  1100 
1725  Jefferson  Davis  Highway 
Arlington,  VA  22202-4136 
Attn:  Dr.  Charles  M.  Huddleston 

Brobeck  and  Associates 

1235  10th  Street 

Berkeley,  CA  94710 

Attn:  Dr.  Francis  C.  Younger 

Chief  of  Naval  Material 
Office  of  Naval  Technology 
MAT-0712,  Roo(ti  503 
800  North  Quincy  Street 
Arlington,  VA  22217 
Attn:  Dr.  Eli  Zimet 

Commander 

Space  and  Naval  Warfare  Systems  Comman 
National  Center  1,  Room  8E08 
Washington,  DC  20363-5100 
Attn:  RADM  Robert  L.  Topping 

Cornell  University 

369  Upson  Hall 

Ithaca,  NY  14853 

Attn:  Prof.  David  Hammer 
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Defense  Advanced  Research  Projects  Agency  C.  S.  Draper  Laboratories 


1400  Wilson  Blvd. 

Arlington,  VA  22209 
Attn:  Dr.  H.  L.  Buchanan 

Dr.  B.  Hui 

Defense  Nuclear  Agency 

Washington,  DC  20305 

Attn:  Dr.  Muhammad  Ova is  (RAAE) 

Dr.  Michael  Frankie 
Dr.  R.  Gullickson 

Department  of  Commerce 
National  Inst. of  Standards  and  Tech. 
Building  245,  B-102 
Washington,  DC  20234 
Attn:  Dr.  Mark  A.  D.  Wilson 

Dr.  Steven  M.  Seltzer 

Department  of  Energy 

Washington,  DC  20545 

Attn:  Dr.  Wilmot  Hess  (ER20:GTN, 

High  Energy  and  Nuclear  Physics) 
Mr.  Gerald  J.  Peters  (G- 256 ) 

Department  of  the  Navy 

Chief  of  Naval  Operations 

The  Pentagon 

Washington, DC  20350 

Attn:  CAPT  T.  L.  Sanders  (0P981N3) 

LCDR  John  Stanovich  (0P981SDI) 
LCDR  Donald  Melick  (0P981SD) 

Dr.  Steve  Bravy  (0P981S0I) 

Mr.  Greg  Montieth 

Directed  Technologies,  Inc. 

4001  Fairfax  Drive,  Suite  775 
Arlington,  VA  22203 
Attn:  Mr.  Ira  F.  Kuhn 

Dr.  Nancy  Chesser 
Dr.  Arthur  Lee 
Ms.  Marla  Shain 

Directed  Technologies,  Inc. 

5945  Pacific  Center  Blvd. 

Suite  510 

San  Diego,  CA  92121 

Attn:  Dr.  Robert  A.  Jacobsen 

Dr.  Harald  0.  Dogliani 

P.  0.  Box  503 

Los  Alamos,  NM  87544 


555  Technology  Square 
Cambridge,  Massachusetts  02139 
Attn:  Dr.  E.  Olsson 

ESL,  Inc. 

Mail  Stop  M-4216 
495  Jova  Drive 
Sunnyvale,  CA  94088 
Attn:  Dr.  Robert  A.  Marth 

FM  Technologies,  Inc. 

10529B  Braddock  Road 
Fairfax,  VA  22032 
Attn:  Dr.  F.  M.  Mako 

GA  Technologies,  Inc. 

P.  0.  Box  85608 
Code  02/503 
San  Diego,  CA  93138 
Attn:  Dr.  Vincent  Chen 
Dr.  Hiroyuki  Ikez 

General  Dynamics  Corporation 
1745  Jefferson  Davis  Highway 
Suite  1000 

Arlington,  VA  22202 
Attn:  Dr.  Daniel  W.  Roth 

General  Dynamics  Corporation 

Pomona  Division 

1675  W.  Mission  Blvd. 

P.  0.  Box  2507 
Pomona,  CA  92769-2507 
Attn:  Dr.  Nen  « .  ridvko 

Mr.  C.  L.  Featherstone 

Grumman  Corporation 
Grumman  Aerospace  Research  Ctr. 
Bethpage,  NY  11714-3580 
Attn:  Dr.  Richard  G.  Madonna 

Headquarters,  Department  of  Army 
DAMOFDE,  Room  2D547 
The  Pentagon 

Washington,  DC  20310-0460 
Attn:  LTC0L  Lou  Goldberg 

HO  Foreign  Technology  Division 
Wright-Pat terson  AFB,  OH  45433 
Attn:  TUTD/Dr.  C.  Joseph  Butler 

HQ  USAF/TXN 

Patrick  Air  Force  Base,  FL  32925 
Attn:  CAPT  Joseph  Nicholas 
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Hudson  Institute 
Center  for  Naval  Analyses 
Alexandria,  V A  22302 
Attn:  Dr.  F.  Bomse 

Hy-Tech  Research  Corp. 

P.  0.  Box  3422  FSS 
Radford,  V A  24143 
Attn:  Dr.  Edward  Yadlowsky 

Idaho  Engineering  National  Lab. 
P.  0.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  Dr.  Francis  Tsang 

Institute  for  Defense  Analyses 
1801  N.  Beauregard  Street 
Alexandria,  VA  22311 
Attn:  Dr.  Deborah  Levin 

Ms.  M.  Smith 

IRT  Corporation 
3030  Callan  Road 
San  Diego,  CA  92121 
Attn:  Dr.  David  Phelps 

JAYCOR 

11011  Torreyana  Road 
P.  0.  Box  85154 
San  Diego,  CA  92138-9259 
Attn:  Dr.  Franklin  S.  Felber 

Dr.  Seung  Kai  Wong 

JAYCOR 

39650  Libery  Street,  Suite  320 
Freemont,  CA  94538 
Attn:  Dr.  Kendal  Casey 

Joint  Institute  for  Laboratory 
Astrophysics 

National  Bureau  of  Standards  and 
University  of  Colorado 
Boulder,  CO  80309 
Attn:  Dr.  Arthur  V.  Phelps 

Kaman  Sciences 
P.  0.  Drawer  QQ 
Santa  Barbara,  CA  93102 
Attn:  Dr.  W.  Hobbs 


Lawrence  Berkeley  Laboratory 
University  of  California 
Berkeley,  CA  94720 
Attn:  Dr.  Edward  P.  Lee 

Dr.  Thomas  Fessenden 
Dr.  William  Fawley 
Dr.  Roger  Bangerter 

Lawrence  Livermore  National  Laboratory 
University  of  California 
Livermore,  California  94550 
Attn:  Mr.  Arthur  G.  Cole 

Dr.  Michael  Delong 
MAJ  Kenneth  Dreyer 
Dr.  Ed  Farley 
Dr.  Alex  Glass 
Dr.  George  Craig 
Dr.  C.  V.  Johnson,  III 
Dr.  George  Kamin 
Dr.  V.  Kelvin  Ne’ 1 
Dr.  Arthur  C.  Paul 
Mr.  Louis  Reginato 
Mr.  Doyle  Rogers 
Dr.  Dennis  R.  Slaughter 
Dr.  David  Whit  turn 
Dr.  Simon  S.  Yu 
Dr.  Frank  Chambers 
Dr.  James  W.-K.  Mark,  L-477 
Dr.  William  Barletta 
Dr.  William  Sharp 
Dr.  John  K.  Boyd 
Dr.  John  Clark 
Dr  Gerrge  J.  Caporaso 
Dr.  Donald  Prosnitz 
Dr.  John  Stewart 
Dr.  Y.  P.  Chong 
Dr.  Hans  Kruger 
Dr.  Thaddeus  J.  Orzechowski 
Dr.  John  T.  Weir 
Dr.  Yu-Jiuan  Chen 

Dr.  James  E.  Leiss 
13013  Chestnut  Oak  Drive 
Gaithersburg,  MD  20878 

Lockheed  Missiles  and  Space  Co. 

3251  Hanover  St. 

Bldg.  205,  Dept  92-20 
Palo  Alto,  CA  94304 
Attn:  Dr.  John  Siambis 


La  Jolla  Institute 
P.  0.  Box  1434 
La  Jolla,  CA  92038 
Attn:  Dr.  K.  Brueckner 
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Los  Alamos  National  Laboratory 
P.0.  Box  1663 
Los  Alamos,  NH  87545 
Attn:  Dr.  L.  Thode 

Dr.  H.  Dogliani,  MS-5000 
Mr.  R.  Carlson,  MS-P940 
Dr.  Carl  Ekdahl,  MS-D410 
Dr.  Joseph  Mack 
Dr.  Melvin  I.  Buchvald 
Dr.  David  C.  Moir 
Dr.  Daniel  S.  Prono 
Dr.  S.  Czuchlewski 
Dr.  Thomas  P.  Starke 
Dr.  Donald  D.  Cobb,  D466 
Dr.  Robert  R.  Karl,  D466 
Dr.  William  B.  Maier 
Dr.  John  P.  Rink 
Dr.  David  Chamberlin 

Maxwell  Laboratories  Inc. 

8888  Balboa  Avenue 
San  Diego,  CA  92123 
Attn:  Dr.  Ken  Whitham 

Dr.  S.  Echouse 

McDonnell  Douglas  Research  Laboratories 
Dept.  223,  Bldg.  33,  Level  45 
Box  516 

St.  Louis,  MO  63166 
Attn:  Dr.  Carl  Leader 

Dr.  Frank  Bieniosek 
Dr.  John  Honig 

Mission  Research  Corporation 
1720  Randolph  Road,  S.E. 

Albuquerque,  NM  87106 
Attn:  Dr.  Thomas  Hughes 

Dr.  Lawrence  Wright 
Dr.  Kenneth  Struve 
Dr.  Michael  Mostrom 
Dr.  Dale  Welch 


National  Inst,  of  Standards  &  Tech 
Gaithersburg,  Maryland  20760 
Attn:  Dr.  Mark  Wilson 

National  Inst,  of  Standards  &  Tech 
Radiation  Physics  Bldg.  Room  C229 
Washington,  DC  20234 
Attn:  Dr.  Wayne  Cassatt 

National  Security  Agency 

4928  College  Avenue 

College  Park,  HD  20740 

Attn:  Dr.  Albert  J.  Leyendecker 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 
Attn:  CAPT  James  Fontana 
Mrs.  Teresita  Finch 
Dr.  Rodney  Buntzen 

Naval  Postgraduate  School 
Physics  Department  (Code  61) 
Monterey,  CA  93940 
Attn:  Prof.  John  R.  Neighbours 

Prof.  Fred  Buskirk 
Prof.  Kai  Woehler 
Prof.  Xavier  Maruyama 

Naval  Surface  Warfare  Center 
Dahlgren,  VA  22448-5000 
Attn:  Dr.  E.  M.  Williams 

Mr.  C.  E.  Gallaher 
Mr.  Lawrence  Luessen 
Ms.  Theresa  Houghton 
Dr.  Ronald  J.  Gripshover 
Dr.  S.  L.  Moran 
Dr.  Edwin  Ball 


Mission  Research  Corporation 
P.  0.  Drawer  719 

Santa  Barbara,  California  93102 
Attn:  Dr.  C.  Longmire 

Dr.  N.  Carron 

Mission  Research  Corporation 
8560  Cinderbed  Road 
Suite  700 

Newington,  VA  22122 
Attn:  Dr.  Khanh  Nguyen 
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Naval  Surface  Warfare  Center 
White  Oak  Laboratory 
Code  R-41 

Silver  Spring,  Maryland  20903-5000 
Attn:  CAPT  R.  P.  Fuscaldo 
Dr.  Thomas  A.  Clare 
CAPT  R.  W.  Moore 
Dr.  Ira  Blatstein 
Mr.  Kenneth  Caudle 
Mr.  Carl  Larson 
Dr.  Robert  DeWitt 
Dr.  Ralph  Schneider 
Dr.  Joel  Miller 
Dr.  Stanley  Stern 
Dr.  Omer  Goktepe 
Dr.  A.  L.  Licht 
Dr.  Joon  Choe 
Mr.  David  Demske 
Dr.  Jag  Sharma 
Mr.  W.  M.  Hinckley 
Dr.  H.  S.  Uhm 
Dr.  R.  Fiorito 
Dr.  R.  Stark 
Dr.  H.  C.  Chen 
Dr.  D.  Rule 
Dr.  Matt  Brown 
Mrs.  Carolyn  Fisher  (G42) 

Dr.  Eugene  E.  Nolting  (H23) 

Naval  Technical  Intelligence  Center 
Code  DA52 

4301  Suitland  Road 
Washing  on,  DC  20395 
Attn:  Mr.  Mark  Chapman 

New  Mexico  State  University 
Research  Center 
Box  RC 

Las  Cruces,  NM  88003-0001 
Attn:  Dr.  Leon  J.  Radziemski 

Northeastern  University 
Dept,  of  Elec.  Engineering 
360  Huntington  Avenue 
Boston,  MA  02115 
Attn:  Dr.  Philip  Serafim 

North  Star  Research  Corp . 
c555  Zuni,  S.  E. 

Albuquerque,  NM  87104 
Attn:  Dr.  Richard  Adler 


Oak  Ridge  National  Laboratory 
Health  &  Safety  Research  Div. 

P.  0.  Box  X 
Oak  Ridge,  TN  37830 
Attn:  Dr.  Rufus  H.  Ritchie 
Dr.  0.  Crawford 

Office  of  the  Chief  of  Naval  Operation 

Strategic  and  Theatre  Nuclear  Warfare 

0P-654E4 

The  Pentagon 

Washington,  DC  20350 

Attn:  Dr.  Yong  S.  Park 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 
Attn:  Dr.  C.  W.  Roberson 

Dr.  F.  Saalfeld 

Office  of  Naval  Research  (2  copies) 
Department  of  the  Navy 
Code  01231C 
Arlington,  VA  22217 

Office  of  Under  Secretary  of  Defense 

Research  and  Engineering 

Room  3E1034 

The  Pentagon 

Washington,  DC  20301 

Attn:  Dr.  John  MacCallum 

0SVR 

P.  0.  Box  1925 
Washington,  DC  20013 
Attn:  Dr.  Jose  F.  Pina 

PhotoMetrics,  Inc. 

4  Arrow  Drive 

Woburn,  MA  01801 

Attn:  Dr.  Irving  Kofsky 

Physics  International,  Inc. 

2700  Merced  Street 
San  Leandro,  CA.  94577 
Attn:  Dr.  E.  Goldman 

Dr.  James  Benford 
Dr.  George  B.  Frazier 
Mr.  Ralph  Genuario 

Princeton  University 
Plasma  Physics  Laboratory 
Princeton,  NJ  08540 
Attn:  Dr.  Francis  Perkins,  Jr. 
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Pulse  Sciences,  Inc. 

600  McCormack.  Street 
San  Leandro,  CA  94577 
Attn:  Dr.  Sidney  Putnam 

Dr.  Vernon  Bailey 
Dr.  M.  Tiefenbach 
Dr.  J.  Edighoffer 
Mr.  James  Fockler 

Pulse  Sciences,  Inc. 

2001  Vilshire  Boulevard 
Suite  600 

Santa  Monica,  CA  90403 
Attn:  Dr.  John  R.  Bayless 

R&D  Associates 
301A  South  Vest  Street 
Alexandria,  VA  22314 
Attn:  Mr.  Ihor  Vitkovitsky 

Dr.  Peter  Turchi 

The  Rand  Corporation 
2100  M  Street,  NV 
Washington,  DC  20037 
Attn:  Dr.  Nikita  Wells 

Mr.  Simon  Kassel 


Sandia  National  Laboratory 
Albuquerque,  NM  87115 
Attn:  Dr.  Collins  Clark 

Dr.  John  Freeman/1241 

Dr.  Charles  Frost 

Dr.  Gerald  N.  Hays 

Dr.  Michael  G.  Mazarakis/1272 

Dr.  John  Wagner/1241 

Dr.  Ron  Lipinski/1274 

Dr.  James  Poukey 

Dr.  Milton  J.  Clauser/1261 

Dr.  Kenneth  R.  Prestvich/1240 

Dr.  Kevin  O'Brien 

Dr.  Isaac  R.  Shokair 

Dr.  J.  Pace  VanDevender/1200 

Dr.  J.  T.  Crow 

Dr.  S.  Shope 

Dr.  B.  N.  Turman 

Dr.  C.  Olson 

Dr.  Richard  Adams 

Dr.  Malcolm  Buttram 

Mr.  Charle:  Crist 

Dr.  Susan  Fisher 

Dr.  John  Keizur 

Dr.  Gordon  T.  Leifeste 

Dr.  Raymond  W.  Lemke 

Dr.  Juan  Ramirez 

Dr.  James  Rice 

Dr.  Michael  Wilson 

Science  Applications  Inti.  Corp. 

2109  Air  Park  Road,  S.  E. 
Albuquerque,  NM  87106 
Attn:  Dr.  R.  Richardson 

Dr.  Michael  D.  Havorth 
Dr.  Alan  J.  Toepfer 

Science  Applications  Inti.  Corp. 

5150  El  Camino  Road 
Los  Altos,  CA  94022 
Attn:  Dr.  R.  R.  Johnston 
Dr.  Leon  Feinstein 
Dr.  Douglas  Keeley 
Dr.  E.  Roland  Parkinson 


Science  Applications  Inti.  Corp. 

1710  Goodridge  Drive 
McLean,  VA  22102 
Attn:  Mr.  tf.  Chadsey 
Dr.  A  Drobot 
Dr.  K.  Papadopoulos 
Dr.  William  tf.  Rienstra 
Dr.  Alfred  Mondelli 
Dr.  D.  Chernin 
Dr.  R.  Tsang 
Dr.  J.  Petillo 
Dr.  G.  Bourianoff 
Ms.  K.  tfilson 

Science  Research  Laboratory,  Inc. 

1600  tfilson  Boulevard 
Suite  1200 
Arlington,  VA  22209 
Attn:  Dr.  Joseph  Mangano 

Dr.  Daniel  Birx 

Commander 

Space  &  Naval  tfarfare  Systems  Command 
PMtf-145 

Washington,  DC  20363-5100 
Attn:  CDR  tf.  Fritchie 
Mr.  D.  Merritt 

Space  Power  Institute 
315  Leach  Science  Center 
Auburn  University 
Auburn,  AL  36845=3501 
Attn:  Prof.  M.  Frank  Rose 

Spectra  Technology 
2755  Northup  Way 
Bellevue,  WA  98004 
Attn:  Dr.  Dennis  Lowenthal 
Dr.  Steve  Baughoura 
Dr.  James  Ewing 

SRI  International 
PSO-15 

Molecular  Physics  Laboratory 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 
Attn:  Dr.  Donald  Eckstrom 

Dr.  Kenneth  R.  Stalder 
Dr.  Roberta  Saxon 
Dr.  Jay  Dickerson 


Strategic  Defense  Initiative  Org. 

SDI0/T/DE0 

The  Pentagon 

Washington,  DC  20301-7100 
Attn:  COL  Thomas  Meyer  (DEtfO) 

LTC  Michael  toole  (DEtfO) 

MAJ  J.  Wills 

Dr.  Dwight  Duston 

LTC  Ed  Pogue 

Dr.  Kevin  Probst 

Dr.  Charles  Sharn 

System  Planning  Corporation 
1500  tfilson  Boulevard,  Room  1213tf 
Arlington,  VA  22209 
Attn:  Mr.  James  T.  Lacatski 

Titan/Spectron,  Inc. 

P.  0.  Box  4399 
Albuquerque,  NM  87196 
Attn:  Dr.  R.  Bruce  Miller 

Dr.  John  Smith 

Titan  Systems,  Inc. 

5910  Pacific  Center  Blvd. 

San  Diego,  CA  92121 
Attn:  Dr.  R.  M.  Dove,  Jr. 

Tetra  Corporation 
4905  Hawkins  Street,  N.  E. 
Albuquerque,  NM  87109-4345 
Attn:  Mr-  tfilliam  Money 

University  of  California 
Physics  Department 
Irvine,  CA  92664 
Attn:  Dr.  Gregory  Benford 

Dr.  Norman  Rostoker 

University  of  California 

San  Diego,  CA  92110 

Attn:  Dr.  Marshall  N.  Rosenbluth 

UCLA 

Physics  Department 
Los  Angeles,  CA  90024 
Attn:  Dr.  F.  Chen 
Dr.  C.  Joshi 
Dr.  J.  Dawson 
Dr.  N.  Luhmann 
Dr.  tf.  Barletta 
Dr.  T.  Katsouleas 
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University  of  Colorado 
Dept,  of  Astrophysical,  Planetary 
&  Atmospheric  Sciences 
Boulder,  CO  80309 
Attn:  Dr.  Scott  Robertson 

University  of  Illinois  at  Chicago 

Dept,  of  Physics 

P.  0.  Box  4348 

Chicago,  IL  60680 

Attn:  Dr.  Charles  K.  Rhodes 

University  of  Maryland 
College  Park,  MD  20742 
Attn:  Dr.  J.  Goldhar 
Dr.  V.  Destler 
Dr.  C.  Striffler 
Dr . Hoon-Jhong  Rhee 

University  of  Michigan 
Dept,  of  Nuclear  Engineering 
Ann  Arbor,  MI  48109 
Attn:  Prof.  Terry  Kammash 
Prof.  R.  Gilgenbach 

University  of  New  Mexico 

Dept,  of  Chem.  5,  Nuclear  Engineering 

Albuquerque,  NM  87131 

Attn:  Prof.  Stanley  Humphries 

Commander 

U.  S.  Army  Laboratory  Command 
2800  Powder  Mill  Road 
Adelphi ,  MD  20783-1145 
Attn:  George  Albrecht  (AMSLC-TP-PL) 

U.  S.  Army  Combined  Army  Center 
ATZL-CAG 

Ft.  Leavenworth,  KS  68027-5000 
Attn:  LTC  Orville  Stokes 

Yale  University 
Mason  Laboratory 
New  Haven,  CN  06520 
Attn:  Dr.  Ira  Bernstein 

Director  of  Research 

U.S.  Naval  Academy 

Annapolis,  MD  21402  (2  copies) 


